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Abstract 

We propose the gauged Thirring model as a natural gauge-invariant general- 
ization of the Thirring model, four-fermion interaction of current-current type. 
In the strong gauge-coupling limit, the gauged Thirring model reduces to the re- 
cently proposed reformulation of the Thirring model as a gauge theory. Especially, 
we pay attention to the effect coming from the kinetic term for the gauge boson 
field, which was originally the auxiliary field without the kinetic term. In 3 -|- 1 
dimensions, we find the nontrivial phase structure for the gauged Thirring model, 
based on the Schwinger-Dyson equation for the fermion propagator as well as the 
gauge-invariant effective potential for the chiral order parameter. Within this ap- 
proximation, we study the renormalization group flows (lines of constant physics) 
and find a signal for nontrivial continuum limit with nonvanishing renormalized 
coupling constant and large anomalous dimension for the gauged Thirring model 
in 3-1-1 dimensions, at least for small number of flavors Nf. Finally we discuss 
the (perturbatively) renormalizable extension of the gauged Thirring model. 
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1 Introduction 



The Nambu-Jona-Lasinio (NJL) model of four-fermion interaction is a suc- 
cessful model to study the chiral dynamics. Actually, it has succeeded to extract 
many characteristic features of the chiral dynamics. However, it is usually thought 
that the NJL model lacks the sound basis as a quantum field theory, namely, it is 
(perturbatively) nonrenormalizable. Therefore, the NJL model is regarded as an 
effective quantum field theory which is nevertheless quite important phenomeno- 
logically. 

Indeed, the NJL model in 3+1 dimensions is non-renormalizable in the usual 
perturbation theory in the coupling constant. Nevertheless, it has been estab- 
lished that a class of the 2+1 dimensional four-fermion model, e.g. the Gross- 
Neveu model, is renormalizable in the different expansion scheme, 1/A^/ expansion 
[0. However, there is no evidence that the 3+1 dimensional NJL model becomes 
renormalizable by choosing some clever expansion scheme. Moreover, it has been 
shown on the lattice that the 3+1 dimensional NJL model is not renormalizable 
even in the non-perturbative sense 0. 

It is generally believed that the non-perturbatively as well as perturbatively 
non-renormalizable model leads to trivial (non-interacting free) theory in the 
continuum limit. Here the continuum limit means the procedure of obtaining the 
well-defined finite continuum field theory from the regularized model by taking the 
limit of removing the cutoffs, e.g. taking the a — > (a: lattice spacing) limit in the 
lattice field theory and taking the A ^ oo limit in the model with an ultraviolet 
cutoff A, etc. Hence, if a theory is nontrivial, it is expected that the theory can be 
renormalizable in a certain renormalization scheme. The authors of ref. |Q use the 
following terminology: they call a theory (non-perturbatively) renormalizable if 
the theory is perturbatively non-renormalizable and nevertheless gives a nontrivial 
continuum limit, rather than simply calling it nontrivial. This is indeed consistent 
with the picture of modern renormalization group (RG) in Wilson's sense 0, 
namely, the nontrivial continuum limit is realized on a nontrivial fixed point 
if the cutoff theory approaches to the fixed point in the bare parameter space 
along the renormalized trajectory as a result of having performed the successive 
renormalization steps. 

Recently, another possibility of renormalizability (in the sense defined above) 
for a class of models with four-fermion interaction was suggested in 0. The 
essential ingredient is the introduction of the gauge interaction into the NJL 
model. If the gauge interaction is incorporated, the model is called the gauged 
NJL model 0. The gauged NJL model was initiated by Bardeen et al. 0] 
motivated from the renormalization group consideration of the strongly coupled 
gauge theory in the sense of Miransky [||, H . It has been found that the gauged 
NJL model has nontrivial phase structures p!0| , |ll| . In the work , it was claimed 
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that the gauge interaction might promote the trivial NJL model to an (non- 
perturbatively renormalizable) interacting theory. It was suggested in [T^ 
that the gauged NJL model is renormalizable if and only if the equivalent gauge- 
Higgs- Yukawa system approaches the infrared fixed point of Pendleton and Ross 
[0 (PR) in the continuum limit and that, for the renormalizability of the gauged 
NJL model, the presence of asymptotically free (or fixed coupling as a standing 
limit of slowly running gauge interaction is essential, but the asymptotic 
freedom of the gauge interaction should not be too strong. Such observations 
were also claimed by Krasnikov in the study of the gauged NJL model in 
4 — e dimensions. This situation is in sharp contrast with the Abelian-gauged 
NJL model with the non-asymptotic free (running) gauge coupling. In such a 
class of the gauged NJL model, there is no indication of the nontrivial fixed point 



on which the nontrivial continuum may be obtained |jT5|, |T6|, |T7|. 



Actually, the renormalizability of the gauged NJL model in 3+1 dimensions 
has been shown within the ladder approximation of the Schwinger-Dyson (SD) 



equation ||18|. And this conclusion has further been supported by the subsequent 
work [jl9| in which the flow of the renormalized Yukawa coupling was analyzed 
in the corresponding gauge- Higgs- Yukawa theory. These works are, however, re- 
stricted to the fixed coupling case in the ladder approximation of the SD equation. 
In the quite recent work of Kyoto group lead by Kugo 0], it was proved in the 
scheme of l/N^ expansion that the gauged NJL model (as a limit of a gauge- 
Higgs- Yukawa system which is perturbatively renormalizable) is renormalizable 
in 3+1 dimensions in the leading order of the 1/Nc expansion. It is desirable to 
extend this analysis to higher orders of l/N^ expansion. 

On the other hand, there is another class of four-fermion interaction, the 
Thirring model pO|, ^ . This model is also important as the low energy effective 



theory of QCD p2|, as well as the field theoretical model for studying the non- 
perturbative aspects of quantum field theory, see references in ||2^, The 
Thirring model is perturbatively non-renormalizable in 3+1 dimensions. How- 
ever, the 2+1 dimensional Thirring model is also renormalizable in 1/A^/ expan- 
sion, as well as the 2+1 dimensional NJL (or Gross- Neveu) model. In light of the 
above argument, it is interesting to pursue the gauge-invariant generalization of 
the Thirring model, which we call the gauged Thirring model. Quite recently, the 
Thirring model was itself reformulated as a gauge theory [^, There 
it has been realized that the gauge structure is actually involved in the original 
Thirring model and the kinetic term for the gauge field which is absent in the orig- 
inal Lagrangian is generated through radiative corrections only when the fermion 
is massive. In this paper we propose a gauged Thirring model. It should be noted 
that the original Thirring model is reproduced as the strong gauge coupling limit 
of the gauged Thirring model, in sharp contrast with the gauged NJL model in 
which the gauge interaction is introduced by hand and the NJL model is obtained 
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as the weak gauge coupling limit. As the gauge- invariant generalization of the 
four-fermion interaction, the gauged Thirring model seems to be more natural 
than the gauged NJL model. In this paper we study the implication coming from 
the gauge-invariant extension of the Thirring model, especially paying attention 
to the effect of the kinetic term for the gauge field. 

We also discuss the possible perturbatively renormalizable extension of the 
gauged Thirring model. It is easy to see that the gauged Thirring model is equiv- 
alent to the gauged non-linear a model with fermions, which implies that the 
gauged Thirring model is perturbatively non-renormalizable. We can extend the 
(perturbatively non-renormalizable) gauged Thirring model into the perturba- 
tively renormalizable model, the generalized gauged Thirring model. However, 
it is not yet clear under what conditions the gauged Thirring models give well- 
defined continuum limits which are equivalent to specific nontrivial generalized 
gauged Thirring models. Such a viewpoint is quite necessary to consider the ex- 
tension of the work 0] to the gauged Thirring model and give an answer for the 



question raised in ref. 



We have found that the recent proposal of gauge-invariant extension of four- 
fermion model by Jersak and his collaborators |2^, |3^ is essentially the same as 
the gauged Thirring model which we have proposed. Their model defined on a 
lattice which they call the xU<P model does not exactly coincides with ours, since 
they use the compact formulation for the gauge field and the staggered fermion 
on the lattice. However, the naive continuum limit of xU<P model reduces to 
the gauged Thirring model, if the staggered fermion is replaced with the Dirac 
fermion. Therefore, their extensive studies of Monte Carlo simulation may shed 
more light on the non-perturbative features of the gauged Thirring model from 
somewhat different directions, although the direct comparison of our result with 
theirs may not be possible. 

This paper is organized as follows. The paper is divided into three parts. In 
part I, we give the formulation of the gauged Thirring model. In part II and 
HI, we analyze the Abelian-gauged Thirring model by making use of the effective 
potential in part II and the Schwinger-Dyson (SD) equation in part HI. 

In part I, section 2 and 3, the gauged Thirring model is formulated. In section 
2, we define the gauged Thirring model so that the Thirring model which is 
recently reformulated as a gauge theory is recovered as the strong gauge-coupling 
limit /? ^ 0. The original Thirring model is regarded as a gauge- fixed version 
of the extended gauge theory, which we call the master gauge theory. Moreover, 
we give the covariantly gauge-fixed BRST formulation of the gauged Thirring 
model. To solve the Schwinger-Dyson (SD) equation in part HI, we need to take 
the nonlocal covariant gauge fixing It is shown that the above formulation 

can be extended to the case of the nonlocal gauge fixing. In sections, we restrict 
ourselves to the Abelian-gauged Thirring model and discuss the limiting cases. 
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Based on the considerations so far, we conjecture the possible phase diagram of 
the gauged Thirring model. In the Abelian gauged case, we can integrate out 
the gauge boson field exactly and can obtain the effective fermionic theory of 
the gauged Thirring model, which reduces to the original Thirring model in the 
strong gauge-coupling limit. This is nothing but the fermionization of the gauged 
Thirring model. However, the effective fermionic theory obtained in such a way 
is nonlocal and has no gauge invariance. 

In part II, from section 4 to 6, we obtain the effective potential for the chiral 
order parameter which is used to obtain the phase diagram of the gauged Thirring 
model. Section 4 is a preparatory section to obtain the effective potential. In 
order to obtain the gauge-independent effective potential, we adopt the inversion 
method invented by Fukuda [03 and developed with his collaborators 133, 



36|| . In section 5, we explicitly obtain the effective potential for the chiral order 
parameter according to the inversion method. It is used to show the existence 
of the chiral phase transition which is the second order. In section 6, from the 
effective potential in the previous section, we obtain the critical line of the chiral 
phase transition in the phase diagram and determine the critical exponent for the 
chiral order parameter. 

In part III, from section 7 to 12, we discuss the phase structure in more detail 
and examine the continuum limit based on the solution of the SD equation for 
the gauged Thirring model. In section 7, we set up the SD equation for the 
gauged Thirring model. We argue the consistency of the ladder approximation 
from the viewpoint of gauge invariance. This is necessary and also possible, since 
the gauged Thirring model has now been formulated as a true gauge theory. 
We show that the bare vertex approximation is consistent with no wavefunction 
renormalization for the fermion. Such a situation is realized only when we take 
the nonlocal gauge. In section 8, we obtain the Cornwall- Jackiw-Tomboulis 
(CJT) effective action whose variation with respect to a full propagator yields the 
SD equation for the full propagator (in the nonlocal gauge). We show that the 
nontrivial (chiral-symmetry breaking) solution for the SD equation gives lower 
effective potential than the trivial (chiral symmetric) solution. This implies that, 
in the chiral limit, the vacuum favors the chiral symmetry breaking rather than 
the chiral symmetric case. 

In section 9, we discuss the solution of the SD equation in the nonlocal gauge. 
We obtain the explicit solution in two regions in the phase diagram which are 
neighborhood of the two critical points of the pure Thirring4 model and the pure 
QED4. In section 10, we write the equi-correlation-length line for the fermion 
based on the solution of SD equation. The critical line is obtained in the limit of 
infinite correlation-length. We discuss also the possible existence of the solution 
in the negative four- fermion coupling region. In section 11, we study the RG flows 
(i.e. lines of constant physics) for the constant renormalized coupling constant 
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and for the constant mass ratio between the fermion mass and the gauge boson 
mass. We examine the continuum hmit of the gauged Thirring model with cutoff 
and find the nontriviality of the continuum gauged Thirring theory. In section 
12, the anomalous dimension for the composite operator -i/^iIj is calculated in 
the gauged Thirring model. It is shown that the anomalous dimension of the 
continuum gauged Thirring theory takes the large value (greater than 1 and 
smaller than 2), if we require the nontriviality of the continuum limit, namely, 
the continuum limit is taken along the RG trajectories obtained in the previous 
section. The final section is devoted to conclusion and discussion. 

In Appendix A, we include the short note for the inversion method which is 
sufficient to give the inversion formula. 
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2 Thirring model as a gauge theory 



We consider the general Thirring model whose Lagrangian is given by 

= ^^^^^^,iJ' -moi:^^^ -^^^Oj^T''^^)^^''^f^T''^p'^)^ (2.1) 

where ip^ is a 4-component Dirac fermion with an flavor index j which runs from 
1 to Nf (summed over from 1 to Nf if repeated) and 7^(/i = 0, 1,2, ...,D — 1) 
are 4x4 gamma matrices satisfying the Clifford algebra {7^,71/} = 2g^i,l, and 
T"(a = l,...,dimG) makes up a basis of the Lie group (color group) G, e.g. 
dimG = N^-l for G = SU{N^). 

By introducing the auxiliary vector field := gA^^T"", the theory is equiva- 
lently rewritten as 

where we have defined the covariant derivative: 

D^[A] := - lA, = d,- igAlT^ (2.3) 
and introduced another parameterization p6| : 



Gr := (2,4) 



The theory given by the Lagrangian ( |2.2| ) is identical with the massive Yang- 
Mills theory where the Yang-Mills field couples minimally to fermions, although 
it lacks the kinetic term for the Yang-Mills field. Note that the kinetic term for 
the fermion and the fermion mass term are gauge invariant, but only the mass 
term for the vector field breaks the gauge invariance explicitly. Therefore, the 
vector field A^ is not still a gauge field at this stage. 

2.1 recovery of gauge invariance 

As is well known from the study of massive Yang-Mills theory the theory can 



be rewritten into the gauge-invariant form by introducing scalar modes 6'" (a = 
1, ...,dimG). Actually, introducing the vector field K^{9): 

K,{e) = -Ued^Ul := ir^(^)T^ Kl{e) = jtT[T'^Ued,Ul], (2.5) 

made from an element of the unitary group U{Nc): 

Ue{x) = exp[i^r(x)T'^] := f/°(x) + tU''{x)T% (2.6) 
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we can write the gauge-invariant generalization of Eq. (pT 



where we have used the following normalization: 

tr(r'^T^) 



(2.7) 



(2i 



The scalar field 6 := g6°'T°' introduced above is nothing but the non-Abelian 
generalization of the Stiickelberg field Q in the Abelian massive gauge theory 
2^]. Indeed, the Lagrangian Eq. ( p.7| ) is invariant under the (finite) local gauge 



transformation: 



U^{x)il)j{x), 

U^{x)A^{x)Ul{x) + -U^{x)d^Ul{x) 



UJx) 



A^{x) - -Ul{x)d^U^{x) 
9 



X) 



-U\x)d^U\x)^ {U\x) := U^{x)Ue{x)) 
U^{x)K,{x)Ul{x) + -U^{x)d,Ul{x), 



UJx) 



K^{x) - -Ul{x)d^U^{x) 



Ulix). 



(2.9) 



/^From the reason explained in the introduction, we further introduce the 
kinetic term for the gauge field which is itself gauge-invariant: 



(2.10) 



where the field strength is defined by := gF^JT°- = d^A^ — d„A^ — ig[A^, A^] 
and \T°-,T^] = ifahcT^ with the structure constant fahc of G. Thus we arrive at a 
gauge-invariant generalization of the Thirring model: 



/5 



2tr(F,.F^'^), 



/5: = 
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2 ■ 



(2.11) 



^This field is also identified with the Batalin-Fradkin (BF) field js^l in the general canonical 
Hamiltonian formalism of Batalin-Fradkin- Vilkovisky |39 (BFV) for the constrained system, as 
pointed out in the previous paper |2^, see the references for explicit derivation of Eq. (2.11) 
from this point of view. 

^ Here we have rescaled the gauge field A° and the scalar field O'^ by g^^ . Hence, in what 
follows, AP should be identified with G^^: AP = G^^. 
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The f3 = case is the Thirring model. We call the model given by the Lagrangian 
Eq. (|2.11|) with non-zero /3 the gauged Thirring model. 

The original Thirring model is identified with a gauge-fixed version of the 
gauge theory with the Lagrangian, Eq. which we call the master gauge 

theory. f\ Indeed, the Lagrangian Eq. ( p.llj ) reduces to the Lagrangian Eq. ( p.2|) , 
if we take the unitary gauge: 9{x) = (and (3 = 0). However, actual calculations 
such as loop calculations are generally impossible in the unitary gauge. For 
such purposes the covariant gauge is most convenient, although both gauge- fixing 
should give the same results on the gauge- invariant quantities, e.g., the chiral 
condensate Various advantages of maintaining such a gauge symmetry 

were emphasized in the previous works 0|. 



The theory with the Lagrangian Crh" can be cast into the form which is similar 
to the Higgs-Kibble model (or gauged non-linear sigma model) in the presence of 
fermions: 

+G-hT[iD,[A]<^yiD''[Am - ^tT{F,,Fn. (2.12) 



with a local constraint: 



$(x)$t(x) = 1, (2.13) 



for the scalar field $(x) := U{x). It should be noted that the theory with this 
Lagrangian is perturbatively non-renormalizable, as well as the original Thirring 
model. 

In the chiral limit mo — *• 0, the gauged Thirring model as well as the original 
Thirring model has the chiral-symmetry, namely, the Lagrangian is invariant 
under the chiral transformation: 

ij^ix) I — y ij'jix) := e'^'^tfjjix), 

^^{x) I — > V^;(x) := V^j(x)e^^^^ (2.14) 



2.2 BRST quantization 

In the following of this paper, we adopt the covariant gauge fixing for the gauged 
Thirring model. In what follows, therefore, we give the covariantly gauge-fixed 

^ The master gauge theory is not unique. Actually we can write another type of the master 
gauge theory for the Thirring model |^^, which is easily extended to the gauged Thirring 
model. This type of the master gauge theory plays the important role of the bosonization of 
the Thirring model. 
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BRST formulation of the gauged Thirring model. It is well known that the total 
Lagrangian 

CtW" = ^Th" + ^GF + ^FP, (2.15) 

which is invariant under the BRST transformation 6b is obtained by adding 
the gauge-fixing term and the Faddeev-Popov (FP) ghost term Cqf+fp to the 
Lagrangian Crh"'- 

CGF+Fp[A,e,C,C,B] = -z55(C'"(F'^[A^]-^5")), (2.16) 

where 5b denotes the BRST transformation. The nilpotent BRST transformation 
5^* = is given by 





= d;\a)c\xI 


SbB''{x) 


= 0, 


6bC%x) 


= -^r"'c\x)c\x), 


SbC%x) 


= -tB\x), 


6biP^{x) 


= z^C"(x)TV(x), 


Sbk;,{x) 


= d;\k)c\x), 



(2.17) 
where 

Df[A]:=6,,d,-zgr'^A';,. (2.18) 
Note that the BRST transformation of can be rewritten into that of 6: 

5Be\x) = K-\eY^c\ (2.19) 

where K-^^Of^ is the inverse of K{9)°'^ which is defined through 

Ki{e) = K''\e)d^e\ (2.20) 

Note that the BRST invariance of the total Lagrangian is guaranteed, since 
Cgf+fp is invariant under the BRST transformation due to nilpotency of BRST 
transformation Eq. (|2.17|) and the original Lagrangian Cth" is also invariant. 
For the BRST transformation Eq. ( |2.17] ), we have 

Cgf+fp = B''F''[A,e] + ^{B^)' - zC''6bF''[A,9]. (2.21) 
After eliminating the B field, we obtain 

^GF = -^{F-[A,e])', (2.22) 
Cfp = -iC''6BF''[A,e] 

= -^C'^ [^^^D;'[A^^^^ (2.23) 
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2.2.1 a covariant gauge 

The conventional covariant gauge-fixing is obtained for tlie following F"". 

F''[A,e]= ^^'Al. (2.24) 

This covariant gauge leads to the following ghost term: 

= -iC^Df[A]C\x). (2.25) 

In this covariant gauge, we notice that there is a crossing term A^K^{6) in the 
Lagrangian Cth"- This term sometimes makes the analysis somewhat cumber- 
some. 



2.2.2 gauge 

We consider another choice for the gauge-fixing function. In the following we 
choose the generalized gauge: 

F^[A, 6] = d^A^ + ^ne), (2.26) 

such that the crossing term A^^K^{6) is canceled with the relevant term in the 
gauge-fixing term: 

Cgf = -y^{F'^[A,e]r = -y^{d^A;y - ne)d^Ai - ^-{rm\ (2.27) 

in the total Lagrangian Cxh'"- For such a cancellation to occur, we must choose 
such that f"' satisfies the following differential equation: 

dT{e) = M^K^ie) = M^K^\e)d^e\ (2.28) 

Instead, the R^ gauge generates a new type of interaction between the ghost field 
and the scalar mode in the FP ghost term: 

Cfp[A, C, C, 0] = -iC^d^Dl\A)C^ - iCX^-^^K-\ey^C\ (2.29) 

In the Abelian case, K""^ = 6ab, we have 

K,{e) = d,e. (2.30) 

Then the equation Eq. ( p.28| ) is easily integrated to give the solution: 

f{e) = M^e. (2.31) 
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This recovers the usual gauge 

F[A, 6] = d^'A^ + iMO. (2.32) 
In the Abehan case, the total Lagrangian is decomposed as ||24[| : 
i^Th'" = ^i>,A + Cg + Cfp, 



^0 = ^id,(^?-kmy 



Cpp = CFp[C,C] = -ic(^d^d>' + -^^C. (2.33) 

In the Abelian case, therefore, the scalar mode can be decoupled for an arbitrary 
^ by choosing the gauge, since the FP ghost term Cpp does not involve the 9 

field, g 

In the non- Abelian case, also, the equation Eq. ( p.28| ) may have the solution, 
since defined by Eq. ( |2.5D satisfies the integrability condition P7| : 

K^^u = d,K^. - d.K + gr^'KlKl = 0, (2.34) 

which is nothing but the pure-gauge condition (or a null-curvature equation). 
Assuming that f°'{9) is expressed in terms of 9 by solving Eq. ( |2.28|) , we can 
write the total Lagrangian: 

C-Th'" = ^H-i''D^[A]i)^ - mQtlj^i)^ 

-^tr{F,^Fn + —{A;r + 5^(9^;:) + ^{B-y 

+—{K;{9)f - \{nO)? + i^FP[A, C, C, 9i (2.35) 

where we have introduced a new Lagrange multiplier field B°'{x) 

In the Non- Abelian case, the scalar mode 9°" is not decoupled, unless we impose 
further condition. Then, in the non- Abelian case, the self-interacting scalar field 
is decoupled only in the very special gauge. More details on the non-Abelian case 
will be given in subsequent papers |^ . 



^ Note that C^.a is aparantly equivalent to the starting Lagrangian of Gomes et al. |f41| . 
However, is a natural consequence of the full gauge-invariant formulation of the (gauged) 
Thirring model and its meaning is completely different from theirs. 
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2.3 nonlocal gauge fixing 

In order to solve the SD equation in the ladder approximation in section 9, we 
need to take the nonlocal gauge. Here the nonlocal gauge is used in the sense 
that the gauge-fixing parameter ^ gets momentum-dependent, i.e. ^ becomes a 
function of the momentum: ^ = C,{k'^)- In the configuration space, the gauge 
fixing term in the nonlocal gauge is given by 

Cgf = —F^lAix), e{x)] J e{y)], (2.36) 



which has been derived for the covariant gauge |^ and the extension to the 
gauge is straightforward, see e.g. ||2^. Here it should be noted that ^"H^^ 
the Fourier transform of ^~^{x), 



IS 



(27r) 



X) 



(2.37) 



while C{k'^) is not the Fourier transform of ^(x), see ref. [^. If ^(A;^) does not have 
the momentum-dependence, i.e., i{k^) — >■ then ^~^{x — y) ^ 5{x — y)^~^ and 
the nonlocal gauge- fixing term reduces to the usual gauge- fixing term, Eq. ( p.22|) . 
In the formal way, the BRST formulation can be easily extended into the case 
of nonlocal gauge-fixing by replacing the product with the convolution when the 
gauge-fixing function ^{x) appears, as shown explicitly for the Abelian case in 



2.4 renormalizable extension 

I want to regard the gauged Thirring model with a special limit of the extended 
theory which is perturbatively renormalizable. The simplest way is to add the 
potential term of the scalar field $ into the total Lagrangian: 

Cpot = /iV($t$) + A(tr(<l>t$))2, (2.38) 



at the same time of removing the constraint Eq. ( p.l3|) . The theory with the 
Lagrangian Crh" + j^pot is equivalent to the Higgs-Kibble model with fermions. 
Therefore, it is quite interesting to study the relationship between the gauged 
Thirring model and the renormalizable extension. This will be done in a subse- 



quent paper 42 
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3 Abelian-gauged Thirring model 

3.1 reformulation as a gauge theory 

In this paper we study in detail the Abehan-gauged Thirring model with the 
following Lagrangian: 

L = ^^{tYd, - mo)^^ - A.^^rr + 7^\id, + tA,)(j)\' - ^F^^F,,, (3.1) 

where j = l,...,Nf and (p is the complex scalar field of unit length |0(a;)| = 1 
parameterized as 

= e*^("). (3.2) 
By using the variable 6, this model is equivalently rewritten into the following 



one: 

C = ^P^{^Y^^ - mo)V'^' - A^^P^Y^^ + ^{A, - d,9f - ^F^^F^.. (3.3) 

This theory is invariant under the local U{1) gauge transformation: 

^Pix) V(a;)e'"^"^ 
A^{x) A^{x) + d^uj{x), 
(f){x) 0(x)e*^(^' {e{x) ^e{x)+iu{x)), (3.4) 

Therefore we need a gauge-fixing term to quantize this theory as discussed in the 
previous section. 

3.2 limiting cases 

In the strong gauge-coupling limit /5 — > 0, this model reduces to the gauge 
invariant-reformulation of the (massive) Thirring model where 6 is nothing but 
the well-known Stiickelberg field. Indeed, after integrating out the bosonic field 
Afj^ under the unitary gauge ^ = 0, we arrive at the Lagrangian of the Thirring 
model: 

£ = i:^{zj>^d^ - mo)^^' - ^(V^^t'^V^'Y- (3.5) 

Therefore the Thirring model is identified with a gauge-fixed version of the master 
gauge theory ([Ol ) or ( p.3|) . From this point of view, it is possible to regard 



the field 6 as the Batalin-Fradkin field in the general theory of the constrained 



Hamiltonian system |39]. In the case of /3 = 0, the field A^ is an auxiliary vector 
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field which is introduced so as to hnearize the four-fermion interaction. Hence it 
does not have the corresponding kinetic term. In our model we have introduced 
the kinetic term, jF^^^F^^, from the beginning, after the identification of the 
field with the gauge field. Therefore, at least for non-zero /3, the field 
is no longer auxiliary but is regarded as true gauge field which obeys the gauge 
transformation, Eq. ( |3.4| ). 

Even in the case of /5 = 0, it has been shown that the kinetic term for the field 
A^j, is generated due to radiative corrections [^, only when the fermion mass is 



dynamically generated in the chiral limit mo = 0. In the case of massive fermion, 
particularly very massive fermion, the generation of the kinetic term can be seen 



also through the procedure of bosonization [^4|, ^ . To see the effect coming 
from the dynamical gauge field A^, we have here introduced the kinetic term from 
the beginning. This corresponds to enlarging the bare parameter space in accord 
with the concept of Wilsonian RG. Our model is regarded as a gauge- invariant 
generalization of the Thirring model where all the fields have the corresponding 
kinetic term. We call this model the (Abelian-) gauged Thirring model. Such a 
proposal, i.e., gauge-invariant generalization of the four-fermion model has been 



given quite recently on the lattice by Jersak et al. |29|. The correspondence of 
this model to ours will be discussed in the final section. 

In the limit 1/Gt 0{Gt — ^ ±C)o), the gauged Thirring model reduces to 
the quantum electrodynamics (QED) with Nf flavors of fermions: 

J^QED = i'K^Yd, - m^W - A,i,^Y^^ - (^F^'F^, (3.6) 

In the weak gauge-coupling limit /? — oo, on the other hand, the gauged 
Thirring model approaches to the non-linear a model: 

Cnl. = ^\d,<l)\' (101 = 1), (3.7) 

since in this limit the gauge field reduces to the pure gauge, A^ = d^'d and 'd can 
be absorbed into the scalar and the fermion field and hence the fermionic part 
reduces to the free theory. 



3.3 order parameter and phase transition 

Here we mention the order parameter which is used to probe the phase transition 
of the gauged Thirring model. In the massless fermion, mo = 0, the chiral conden- 
sation {ipip) can be used to distinguish the spontaneous-chiral-symmetry-breaking 
phase from the chiral-symmetric phase. For the chiral-symmetry-breaking phase 
transition, we know that there exists one (non-zero and finite) critical coupling 
in the Thirring model as well as the coupling constant j3c in QED4. Therefore it 
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will be natural to imagine that there will exist a critical line which connects two 
critical points and the chiral symmetry is broken spontaneously below this critical 
line in the phase diagram (/?, Gt)- The expected phase diagram is schematically 
shown in Fig. 1. In this paper we will confirm this picture and obtain the explicit 
form of this critical line based on the effective potential for the order parameter, 
(ipip) as well as the solution of the SD equation for the fermion propagator. At 
the same time, we determine the order of the chiral phase transition together 
with the critical exponent. 

3.4 gauge fixing and fermionic action 

In the following we discuss the behavior of the Lorentz-covariant and gauge- 
invariant order parameter, especially the chiral condensate. Since it is gauge 
invariant, it should not depend on which type of covariant gauge-fixing we may 
adopt in the quantization. In this paper we choose the most convenient gauge 
in performing the calculation. In the previous section, we have shown that, if 
we choose the gauge, the Stiickelberg field 9 in the gauged Thirring model is 
completely decoupled irrespective of the value of the gauge-fixing parameter ^. 
Then the total Lagrangian Ctot is decomposed into C^^a and Cg, see Eq. ( |2.33[ ). 
The Lagrangian C^^a is rewritten as 

+Ia^{x)D^^J-\x)A^{x), (3.8) 

where 

D^'t\^) = i-f3d^d, + M^)g,, + {P + r')d,d,- (3.9) 

This is the inverse of the bare gauge boson propagator in the gauge. In 
momentum space it reads 

= f^k' (^g,. - ^) + M'g,^ - r%K. (3.10) 

In the Abelian gauge theory, we can completely integrate out the gauge field 
Afj^. Then we get the effective fermionic action f\oi the Abelian- gauged Thirring 
model: 

= ^P'ix)itrd, - mo)^\x) + ^ip\x)r^^ix)D^^X^ " y)^p\y)Y^p\y),i^.n) 

^The summation over the spacetime point is understood if it is repeated. 
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where the gauge boson propagator is obtained as 



In the strong gauge-coup hng hmit /5 = 0, we have 

(3.13) 

Then it turns out that the gauged Thirring model in the gauge reproduces the 
original Thirring model in the limit ^ ^ oo (namely, the unitary gauge ^ — >■ 0), 
since 

DfKx -y)^ M-^g.J^'ix - y) = Grg^j'^ix - y). (3.14) 

On the other hand, in the weak gauge coupling limit /5 +oo, S'o reduces to 
the action of the free fermion. Note that the effective fermionic theory given by 
Eq. (|3.11| ) is nonlocal and has no gauge invariance, while the chiral symmetry is 
preserved in the chiral limit. 



3.5 renormalizable extension 

It is interesting to consider more general model where the scalar field is not fixed 
to be of a fixed length: 

+ - ^"M" - ^d'^l')'' (3-15) 

Separating the scalar field as follows, 

0(x) = -l=p{x)e''(^\ (3.16) 

we have 

\{d,-tA,m' = ^{d,pY + ^p'{d,9r + ^p'A,A^ - p'A.d^e. (3.17) 

^If we require the current conservation: d^{'4j{x)^^^{x)) — 0, the fermionic action does not 
depend on the gauge-fixing parameter. 
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In the unitary gauge 6 = 0, this term reads 

\{d, - zA,)<p\' = lid.pf + ^p'A.A^, (3.18) 

If we fix the length of the scalar field: p{x) = po, then we have 

\{d, - iA,)<p\' = ^pl{d,ef + ^plA,A^ - plA,d^e, (3.19) 

This is a renormalizable extension of the gauged Thirring model. Therefore, it 
is interesting to perform the renormalization group (RG) analysis of the model 
and study the relationship with this model and the gauged Thirring model. This 
will enable us to extend the phase structure of the gauged Thirring model and 
consider the universality between this model and the Thirring model based on 
the gauge-invariant formulation. The gauged Thirring model is obtained as the 
limit A ^ oo from this model by adjusting the bare mass parameter ap- 
propriately. Hence the usual perturbation theory with respect to the coupling 
constant A is insufficient to study the gauged Thirring model and we need some 
non-perturbative methods. This issue will be discussed in subsequent papers. 
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4 Chiral order parameter 



In the following three sections, we pay attention to the chiral order parameter 



1 



N 



/ 7 = 1 



(4.1) 



and obtain finally the effective potential in terms of the chiral order parameter 



according to the inversion method 



see Appendix. This order pa- 



rameter is gauge-invariant. Therefore this quantity should not depend on which 
type of gauge fixing we may choose. In the following, we choose the gauge 
which is the simplest gauge to calculate the chiral order parameter in the gauged 
Thirring model, since in this gauge the scalar field is completely decoupled irre- 
spective of the gauge-fixing parameter. 

In order to evaluate the chiral order parameter and discuss the spontaneous 
chiral symmetry breaking in the chiral limit tuq = 0, we introduce the follow- 
ing local source term into the Lagrangian according to the rule of the inversion 
method: 



(4.2) 



It is known that the vacuum expectation value of the fermion composite op- 
erator %l){x)ip{x) is calculated according to the expression: 



d"k tr(l) J ^ _ ^1 



(27r)^p2 + J2 2N^ 



f 



(2vr) 



D 



dJ 



(4.3) 



if we use the following form of the gauge boson propagator: 



(4.4) 



This relation shows that the chiral condensation can be calculated irrespective of 
the longitudinal part of the gauge boson propagator. This implies that the chiral 
order parameter obtained according to Eq. ([4.3|) does not depend on the explicit 
form of Dl and hence the gauge-fixing parameter ^. This is expected from the 
gauge invariance of the chiral condensate. However, in the usual calculation, a 
specific gauge is chosen for calculating the gauge-invariant quantity, see e.g. ^ 



The origin of this difference has been explained in ref. |56]. This form has been 



already used to analyze the gauged NJL and the gauged Yukawa models |47 



the Thirring model in D (4 > D > 2) dimensions ||25| and the three-dimensional 
QED [|^, in order to obtain the gauge-invariant phase structure. 
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The first part of tlie riglit-liand-side (RHS) of Eq. ( |4.3| ) is obtained from the 
following relation for the fermion propagator S: 

, 1 In , In r^n f ^ tr(l)J 

-lndet[S„-] = -t.ln[S„-'l = t.^ = / (4.5) 
which comes from the part: 

lny"[rf^][rf^]exp J d^xtpSo^tfj =lndet5'o"\ (4.6) 
where 5*0 is the bare fermion propagator in the presence of source term: 

So{p) = {l'p, + J)-'. (4.7) 
The second part of the RHS of Eq. ( [4.3|) has an alternative expression: 

/ (Ss^-**'!?"''"**' = / - ^)'='Or(k}§jn(k). (4.8) 

Here we have used the fact that any gauge-invariant regularization leads to the 
transverse form for the vacuum polarization tensor of the gauge boson field: 

n^.(fc) = {9^.uk'' - k^KMk). (4.9) 

The formulation of Thirring model as a gauge theory justifies the use of gauge- 
invariant regularization and guarantees the regularization-independence of the 
result. 

The above expression is slightly different from the formula for the chiral order 
parameter obtained by Ukita, Komachiya and Fukuda |34|. They have taken into 
account up to the lowest order of the series expansion in the coupling constant. 
In the gauged Thirring model the corresponding term in the lowest order is given 

byO 

d^k _/n\/,x d ^ r d^k ,^ ,x,9_m-).,x d 



(2vr) 



nDtKk)^^^'^{k) = J ^{D - l)eDP{k)^U{k), (4.10) 



where 



^ The inversion up to the lowest order in the coupUng constant corresponds to the quenched 
approximation in the following sense. If we adopt the fermion propagator as the nonlocal order 
parameter, the lowest-order inversion method leads to the the SD equation in the quenched 
ladder approximation where the photon propagator is replaced with the bare one as well as the 
vertex function, see Appendix. 
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To obtain the explicit form of the vacuum polarization function Il{k), we here 
adopt the dimensional regularization as a gauge-invariant regularization. This 
leads to the following well-known result: 



-Nf 



2tr(l) 
(47r)«/' 



re'r(2 - - 



da 



a{l — a) 



[J2- A;2a(l-a)]2-«/2' 



(4.12) 



da- 



2 J ail 



a 



For D = 4, we get in the space-like region A;2 = —k"^ < 



U{k) = N 



27r2 



1/ A2 



5\ 



(4.13) 



(4.14) 



(9 p' 



J 



2 /3 72 



kl (P^2 



(4.15) 



see 



48| , ^ for details. The derivative Eq. (|4.15| ) is unchanged even if we take 



other regularizations, e.g, the Pauli-Villars regularization, see although the 
finite constant in n(/c) depends on the regularization scheme adopted. 
Now we move to the Euclidean space. In the lowest order, we have 

(27r)4 ^""^^0.7 ^ ^ 



2Nf 
3 

32^ JO 
4^8^ 



dx 



X 



f3x + M'^ 

1 1 M2 



P /52 A2 



2 73 72 
J + ^ln^ + 0(j5) 

, /3A2 + M2\ ^ 
In — J 



J3 Inf 
+^A^io /3a; + M2 



M2 
C(J^) 



(4.16) 



where we have defined x := k\ 



— A;2 > and introduced the ultraviolet (UV) 
cutoff A. 

Then the chiral order parameter is obtained as a function of the external 
source (fermion mass): 



ip := 47r' 



A3 



1 + 



87r2 



I) P 



K 



J 

A 



3 J3 /-A' ^ lnj2/a; 



(4.17) 
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where we have defined the dimensionless couphng constant k by 

_ M2 1 



(4.18) 



Note that the above effective potential has no dependence on the flavor number 
Nf in the lowest order |BD . 



5 Effective potential for chiral order parameter 



Inverting Eq. ( [4.17| ) with respect to J according to the inversion formula 
see Appendix), we obtain 



J 3 



dt 

Jo 



pt + K 



where we have deflned 



T = TiP,K) :=1 



3 1 



(5.1) 



(5.2) 



In order to obtain the effective potential, we use the following relation (see Ap- 
pendix A) between the effective potential and the translation-invariant source 
J: 



dVjif) ^ J_ 
dif ^ A' 



where we have defined the dimensionless effective potential V{ip) by 

(p := 47r^ 



First, we consider the limit (3—^0. Then we have 



A 



where 



r O,K =1 , := — -. 

K loTT^ 



This reproduces the effective potential of the Thirring model: 



T 2 4 /i 2 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



(5.7) 
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Next, we consider the case of k — 0. This case corresponds to the pure QED4. 
Eq. (|5.1| ) reproduces the result obtained in ref. [0. | 



J 
A 



(5.8) 



where 



r(AO) = l- 



/5 



(5.9) 



87r2 



(5.10) 



Then the effective potential for QED4 reads 



For non-zero (3 and k, we have 
dip " A 



r 2^/3' 
2^ 



(5.11) 



In M + ^ j ln(^2 - PolyLog 



2,- 



/9' 



(5.12) 



where we have introduced the polylogarithm function: 

rO log(l-t) 



PolyLog[2,z] 



t 



-dt=Y: 



00 



2 ■ 



n=l 



n 



(5.13) 



Therefore, the effective potential for the chiral order parameter in the gauged 
Thirring model is obtained for arbitrary j3 and k: 



Iny?^ + 



PolyLog 



.(5.14) 



This is an extension of the work ||3^ for QED4. It should be noted that the k = 
case is somewhat special in the sense that the (y9^(ln yj^)^ term does not appear if 
the QED4 limit is approached from non-zero k, see Eq. ( |5.14|) . In fact, the k = 
case has an ambiguity coming from the lower bound of the integral Eq. ( |5.1| ). 

The order parameter ip can have a nontrivial minimum for r < 0, as long as 
the second p'^ part takes the positive value and convex in p when r = 0. This is 
satisfied for the effective potentials Eq. (|5.11| ) and Eq. (IsTTI) at least in the small 



A 2 

Here we have used dtjln^ = —5 (In 



V \2 
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field region <^ 1, although Eq. ( |5.7| ) is negative for the large field > 0{1), 
Fig. 2 (a), (b). 

The term (p"^ (^—\n(p'^ + is positive for ip <^ 1. However, PolyLog[2, —t] is 
negative for t := jS/n > and the absolute value is monotone increasing in t. 
Nevertheless, there is a certain value ipo > such that the second term in (|5.14|) 
is positive for < < ipo, no matter how large the ratio P/n may be, see Fig. 2 
(c). 

Note that the effective potential ( 5.14 ) obtained in the lowest order does not 
have any Nf dependence. In the original work [Q, the effective potential (|5.11 



has been used to estimate the critical coupling of quenched QED4. However, it 
is well known that the effective potential (|5^ ) is justified in the limit Nf 00, 
see e.g. 0. Moreover, the critical point P'^ of QED4 obtained from the effective 
potential ( 5.11 ) agrees with the result for Nf > 1 obtained from the SD equation 
[Q. In view of these, therefore, it is reasonable to consider that our result ( |5.14| ) 
will be applied to the case of relatively large Nf in the gauged Thirring model. 



6 Phase diagram and critical behavior 
6.1 critical line 

We can specify the location of the critical line by putting r(/5, Gt) = which sepa- 
rates the spontaneous chiral-symmetry-breaking phase from the chiral-symmetric 
phase. If we define the reduced coupling constants: 

/3:=^, k:=-, (6.1) 
yc f^c 



by using the critical couplings (3c of pure QED and Kc of Thirring model, the 
equation of the critical line is written as 

^+l|m + 2f) = l. (6,2) 

In the phase diagram [jS, k) the chiral symmetry is spontaneously broken below 
this critical line. The critical line Eq. ( |6.2|) is plotted in Fig. 3. It is easy to show 
that the critical hne reduces to 

= 1-^/3 + W), (6.3) 

in the neighborhood of (0, k^). | 

^ This shows good agreement with the resuU of Monte Carlo simulation of the corresponding 
model of compact version on the lattice by Jersak et al. p9[. 
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The critical line ( p.2| ) obtained from the chiral order parameter is gauge in- 
dependent, although we adopted the special gauge to simplify the calculation. 
Of course, it should be noted that the other critical lines can exist, if we use 
other order parameters (local or nonlocal). Hence the critical line ( |6.2D is not the 
only one for this model. Actually, by taking into account the order parameter 



containing the scalar field, more rich phase structure will be probed, see |29 



6.2 critical exponent and order of the transition 

The effective potential ( [5. 141 ) shows that the chiral order parameter exhibits the 
power law behavior in the neighborhood of the critical line, apart from the loga- 
rithmic correction: 

(t/'V') ~ kl^'11og|r|P, (6.4) 

where r denotes the deviation from the critical line in the two-parameter space 
{i3,k). Our result is consistent with the mean- field theory and the critical ex- 
ponent Pch for the chiral order parameter takes the mean- field value Pch = 1/2, 
irrespective of the values of the coupling constant. The effective potential for the 
chiral order parameter Eq. ( ^.14| ) shows that the chiral phase transition is of the 
second order on the whole critical line. A possible explanation of this disagree- 
ment with the lattice result for the compact version p9[ will be given in the final 
section. 
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7 Set up of the Schwinger-Dyson equation 



In the following sections, we study the phase structure in detail by the use of 
the Schwinger-Dyson (SD) equation. The effective potential for the chiral order 
parameter obtained in the previous section is gauge independent. However, as is 
well known, the solution of the SD equation for the fermion propagator is never 
gauge independent and depends on the chosen gauge. However, the criterion (the 
renormalized coupling constant or the mass ratio) which we use in section 11 for 
the nontriviality of the gauged Thirring model is gauge independent and does 
not depend on the specific gauge. Therefore, we can hope that the final result 
obtained from the SD equation on the gauge-invariant quantities does not depend 
on the choice of gauge which is necessary to actually solve the SD equation. 



7.1 introduction of the nonlocal gauge 

We can write down the SD equation for the gauged Thirring model. The full 
fermion propagator is written as 

S{p) = [A{p')Yp,-Bip')]-\ (7.1) 

in accord with the bare fermion propagator 

So{p) = {Yp,-mo)-\ (7.2) 



For the gauged Thirring model with the Lagrangian ( |2.33| ), the exact SD equation 
for the full fermion propagator is given by 

S-\p) = S,\p) + J £^j^S{q)r^{p,q)D,^{p-q), (7.3) 

where Ti^{p, q) is the full vertex function and D^y{j) — q) is the full gauge boson 
propagator. Note that this SD equation can be decomposed into a pair of inte- 
gral equations for the wavefunction renormalization function A{p'^) and the mass 
function BijP'). 

The SD equation for the full fermion propagator constitutes the closed set of 
equations together with the SD equations for the full vertex function and the full 
gauge boson propagator. In order to solve the SD equation for the fermion propa- 
gator, therefore, we must anyway specify the full vertex function Ty{p, q) and the 
full gauge boson propagator D^y{p — q). In most cases, an ansatz for the vertex 
function is adopted instead of solving the SD equation for the vertex function, 
although it is in principle possible to solve all the SD equation simultaneously. 
The simplest choice is the bare vertex approximation: 

r^(p,g) = 7M- (7.4) 
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This leads to the simplest ladder approximation. For this approximation to be 
consistent with the Ward-Takahashi identity which is a consequence of the gauge 
invariance of the theory, Q there should be no wavefunction renormalization for 
the fermion. Therefore, the bare vertex approximation should self-consistently 
yield the result: 

A^p") ^ 1, (7.5) 

as the solution of SD equation Eq. ([7I3|). This requirement is satisfied if we in- 
troduce the momentum-dependent function //(/c^) in the gauge boson propagator 
as 



(7.6) 



As is shown in the next subsection, the function r]{k'^) can be chosen so that 
the SD equation Eq. ([7.3|) in the bare vertex approximation leads to the solution 



Eq. (|7.5| ) for the wavefunction renormalization A{p ) |3^. Therefore, in this 
setup, we have only to solve the SD equation for the fermion mass function B{p^) 
in the nonlocal gauge. 

The introduction of momentum-dependent function ?7(/c^) in the full gauge 
boson propagator corresponds to introducing the nonlocal gauge fixing term into 



the total Lagrangian instead of the usual one, see Eq. ( |2.36| ). By comparing (|3.12| ) 



with (|7.6|) , the correspondence between ^ and t] is given in momentum space as 
follows. 

' r]{k^)M^ + pk^ -U{k)' ^ k^ + ^{k^)M^ ■ ^ • ^ 

7.2 SD equation in the nonlocal gauge 

The SD equation is decomposed into a pair of integral equations: 

^ , tr[S(p)7'^p^] 



p^tr(l) 



where S is the self-energy part: 

S(P) := / ^l,S{qh,D,,{p - q). (7.9) 

In order for the vertex function to be consistent with the Ward-Takahashi identity, the 
ansatz for the vertex function should include S, i.e. A and B. However, this requirement is not 
sufficient to determine the vertex uniquely, see ref. |5]| and references therein. 
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The SD equation for the fermion wave function renormahzation A is given by 



fc2 ' \ k'^ k^ 

On the other hand, the SD equation for the fermion mass function B reads 



. (7.10) 



mo + e 



-DT{k)[D-r^{k% {7.11] 



{2tt)D q^A^iq^) + B^iq^y 
Separating the angle {} defined by 

k'^ := {q — p)"^ = X + y — 2^yxy cosi!}, X := p"^, y ■= q"^, ('''•12) 

we find 



xA{x) — X 



f 

Cdc^ / dy 
Jo 



Jo 



xk' 



yA^iy) + B^iy) 
"xy cos i)[D — 2 + rjlk"^)] ^xy — {^fxy cosd) 



fc2 



and 



Bix) = mo + CdS^ / dy 
Jo 



A2 y(,D-2)/2Q^y^ 



yA\y) + B\y) 



k^ 



K{x,y), 



-^(^') , (7.13) 



K{x,y) := r d'dsm'^-^ ^DT{k)\D - 7]{k% 
Jo 



where 



2Dt^{D+1)/2Y(^ 



D-W 



(7.14) 



(7.15) 



(7.16) 



Following the same procedure as given in the Appendix of ref.[^, it turns 
out that the requirement A{p'^) = 1 is achieved if ri{k'^) satisfies the following 
equation: 



z''-^Dt{z)7]{z))' + {D- 2)z''-^D't{z) = 0, z:= 



D-l n' 



(7.17) 
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This is simply solved as follows. 



D-2 



t{z} Jo 



dtD'T{t)t 



(7.18) 



where we have assumed that [z'^'^ DT{z)rj{z)]\z=Q = so as to eliminate the 
l/z^~^ singularity in ri{z). This should be checked after having obtained the 
function ri{z). Alternatively, we can write 



7^{z) = {D-2) 



D - 1 



zD-^Dt(z) Jo 



dWT{t)t 



(7.19) 



where we have assumed that [z D'];'{z)\\z=o = 0. 

Thus, once the function Dxik^) is given, we can find the nonlocal gauge r]{k'^) 
according to the above formulae, Eq. ( |7.18|) or Eq. ( [7.19| ), so that A{k'^) = 1 
follows in agreement with the bare vertex approximation. Then we have only to 
solve Eq. ( |7.11| ) for the fermion mass function BijP'). 



mo + e 



-DT{k)[D-r^{k% 



Note that, if we use the following definition of gauge boson propagator: 



kfj_ky 



DL{k) 



kfi ky 

1^'' 



the integrand is replaced as follows. 

DT{k)[D - r/(P)] = {D- l)DT{k) + DL{k). 



(7.20) 



(7.21) 



(7.22) 



8 CJT effective action in the nonlocal gauge 



8.1 CJT effective potential 



Now we derive the Cornwall- Jackiw-Tomboulis (CJT) action such that the 
stationary point of the CJT effective action as a functional of the fermion prop- 
agator gives the SD equation for the fermion propagator in the nonlocal gauge. 
In the case of QED3, it has been already given in ref. [Q. 
The CJT effective action |53] is given by 



T[S] 



:= -Tr[ln(S-^So) + S^^S - 1], 
1. 



'-Ti[-i^'SYSD 



(8.1) 
(8.2) 

(8.3) 
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where we have chosen Tq[S] so that Tq[S = So] = 0. Indeed, this leads to the SD 
equation at the stationary point (see Fig. 4): 







6T[S] 



S 



-1 



n-l 



Exphcit evaluations in momentum space lead to 



Tr[ln(5-^5o)] tr(l) f d^p 



(2vr) 



D 



In 



1 + 



5V 

p2 



.5) 



nSo'S] ^^^^^^ f d^p B{p')[mo-B{p')] 



n 



(2vr; 



D 



p2 + B2^p2^ 



(8.6) 



where Q is the volume of spacetime, ^ = J d^x and tr(l) denotes the trace over 
the spinor indices, e.g. tr(l) = 4 for D = 4. The second part is also evaluated as 



Ij[S] 



dp f d q 



■Tr[YS{p)YS{q)]D,,{p-q), 



2 J (27r)^7 {27t)D 
tr(l) f d^p f d^q [D -r]{e)]DT{k)B{p'^)B{q'^) 



{2n)D J (27r)^ [p^ + B^{p^)][q^ + B^{q^ 



1.7) 



where we have used the observation that the A-dependent term appearing in 
the integrand vanishes if we adopt the nonlocal gauge obtained in the previous 
section, see Eq. ([7.8|). 

For the translation invariant S, i.e., S{x,y) = S{x — y), we can define the 
CJT effective potential V as 



V[B] 



m 



V[B] 



tiil] 



d^'p 
(2^ 



In 



1 + 



5V 

p2 



2B{p')[B{p') ~ mo] 
p^ + B^{p^) 



tr(l) f d^p f d^q DT{k)[D - 7]{k'^)]B{p^)B{q^ 



(27r)^y (27r)^ [p^ + B^p'Y^lq^ + B^q')] 



If we separate the integration variable into the radial part x := p^ and the 
angular part ■(9, and introduce the UV cutoff for the radial part, the effective 
potential reads 



V[B] = 

2 JO 



1 + 



BHx) 



X 



2B{x)[B{x) - mo] 
X + B'^(x) 



^CoCl iy (8.10) 
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where K{x, y) is the kernel defined by Eq. ( [7.15| ) and we have defined 

1 



C 



D 



c. 



D- 



2^7r^/2r(f^ 



(8.11) 



Indeed, the variation of the CJT effective potential with respect to B, i.e 

6 



-V[B] = 



8.12) 



yields the SD equation Eq. ( |7.20| ). However, it should be noted that this is 
possible only when the functions Dt and rj do not include any 5-dependence. 
Hence the following argument is restricted to the quenched case: H = or the 
restricted (approximated) unquenched case ||15|, |l^ where the massless fermion 
is used in evaluating the vacuum polarization function. 

If we denote the solution of the SD equation by B^oi, the effective potential 
at the stationary point reads 



V[B. 



trfr 



soil 



D 



In 



1 + 



Bsoi{p^)[Bsoi{p^) -mp] 

p' + bUp^) 



.(8.13) 



Note that the SD equation has always a trivial solution when mo = 0. This ex- 
pression shows that the nontrivial solution Bgoi 7^ gives lower effective potential 
than the trivial one Bgoi = 0, i.e., V[i?so/] < V[Bsoi = 0] = 0. This is shown by 
rewriting Eq. ( |8.13| ) as 



V[B^ 



sol\ 



tr(l)^ 



9{t) 



2 

ln(l + t) 



A2 



2/ 2n(D-2)/2 



dp (p 



9 



p2 



t 



1 + t' 



(8.14) 



since ^'(0) = and g{t) is positive and monotone increasing function of t. There- 
fore, if the SD equation has a nontrivial solution in the chiral limit mo = 0, the 
chiral-symmetry-breaking solution should be realized. 



8.2 scale invar iance 



By introducing the scaling parameter /x, we can rewrite the effective potential 
( CT as 



V[B] 



+ /i 



-2D+2 



c 



.2A2 .u2A2 



Blix) 



D 



JO JO \ 



_ 2B^{x)[B^{x) - finio] \ 
x + Blix) J 

x^Bf,{x)y^B^{y) 
V^2'^2y + Bl{x)][y + Bliy)] 



(8.15) 



31 



where we have defined the scaled mass function B^ijp') by 



(8.16) 



We denote the effective potential ( p.lO|) by V[B; A,mo; {««}] where {ai} denotes 
the set of all the other bare parameters appearing in the kernel K{x,y; {ai})- 
Then the effective potential ( |8.1CI| ) satisfies 



V[B; A, mo; {ai}] = /i [B^; fiA, /imo; {a'^}] 



(8.17) 



or equivalently 



y[5^;A,mo;{a:}]=/i^V 



^ A mo 

fj, fj, 



lit 



if the kernel satisfies the following property 



K 



'p2 q2 



(8.19) 



under the appropriate transformation between the new bare parameters a'^ and 
the original bare parameters a^, i.e. a- = fi{{aj}j, fi). 

If we require ( |8.19| ) for the kernel, the scaled mass function satsifies the same 
type of the SD equation as Eq. (|7. 14|) for B{p'^) with the bare parameters mo and 
{ai}: 



B^{x) = fimo + Cd dy 
Jo 



K{x,y;{ai}), 



y + Bliy) 



(8.20) 



but in the scaled range of momentum p"^ G [0, (/iA)^] and with the transformed 
bare parameters, /xmo and {a^}. 

In the quenched QED4, the relation (|8.19|) is trivially satisfied without any 
change of the coupling constant, since the kernel has the form: 



K{x,y) = K{x,y;e) oc 



(x-y) ^ 9{y - x) 



X 



y 



(8.21) 



In the gauged N JL4 model with the fixed gauge coupling []TB| , the kernel is given 
by 



K{x,y) = K{x,y;e,G) oc 



9{x-y) 6{y-x) 



X 



G. 



(8.22) 
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Hence the relation ( p.l9| ) is satisfied when 



(8.23) 



i.e. the gauge couphng constant is unchanged and only the four-fermion coupling 
constant G is scaled. 

For the gauged Thirring model, such a transformation between the bare pa- 
rameters is given by 



(8.24) 



as can be seen in the next section, although the explicit kernel is different from 
the gauged NJL model. However, it should be noted that this does not give the 
exact RG flow in the gauged Thirring model, see section 11. 

When all the bare parameters are unchanged, differentiating both sides of 



Eq. (|8.18|) with respect to /i and putting /i = 1, we obtain 



V[BsouKmQ] = — 



2A' 



d 



+ mo 



d 



dmo 



V[B;A,mo] 



B=B,. 



where we have used 



5_V_ 
5B 



B=B., 



0. 



In applying this formula to Eq. (|8.14|) , note that 

d d 



2K' 



dmo 



B=B, 



Hence, it is wrong that 



2A^ 



d 



+ mo 



d 



dmo 



V [B; A, mo] 
V [B sol] A, mo] . 



V[B,,f,A,mo = 0] = -tr(l)-^A^(7 



CD,njB\A' 



A2 



(8.25) 



(8.26) 



(8.27) 



(8.28) 



r?(t) = ln(l + t) - 



1 + t 



2t^ 3t* 



(8.29) 
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9 Solution of SD equation 

In the following sections we adopt the quenched (ladder) approximation in which 
the vacuum polarization function 11 is neglected in Dt'- 

Dr(k^) - «?'(*^) - (8.1) 

This approximation corresponds to the Nf — > limit. Hence the following result 
will hold for the relatively small Nf, in sharp contrast with the result of the 
effective potential. In the quenched limit, the nonlocal gauge is given by 



ri{k 



2\ 



6 f"' •-(0),.u2 



/ M2 



(A;2)3 Jo (3t + M2 

(9.2) 



1-2— + 2— -In 1 + ^ 



First we observe two limits: In the pure QED limit: M — > 0, rj^k'^) approaches 
a constant uniformly in k"^: 

Vik') ^ 1, (9.3) 

which corresponds to the Landau gauge, ^ = 0. 

In the pure Thirring-model limit: f3 ^ 0, on the other hand, rjlk"^) vanishes: 

r/(fc2) -> 0. (9.4) 

This corresponds to the unitary gauge, ^ ^ oo. The nonlocal gauge is plotted as 
a function of jSk"^ /M'^ in Fig. 5. The function rjik'^) increases monotonically and 
has values between and 1. 

In the nonlocal gauge (|9.2D, we study the solution of the SD equation: 

B{x) = mo + dy ^f^f K{x, y), (9.5) 

lovr^ Jo y + B^[y) 

where the kernel K{x, y) is given by 

K{x,y):=- r dds\Yi^'dDT{k)[A-'q{k% (9.6) 

TT JO 
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and we have defined the squared Euchdean momenta: 



X 



V 



> 0, 



y:=q 



> 0. 



(9.7) 



Here we have introduced the UV cutoff A as an upper bound of integration. 

In the chiral hmit mo = 0, the equation Eq. (|9.5| ) has a trivial solution 
Bip^) = irrespective of values of the coupling constants. If the SD equation 
has a nontrivial solution in the chiral limit mo = 0, then the fermion mass is 
dynamically generated and the chiral symmetry is spontaneously broken, since 
the nontrivial solution gives lower effective potential than the trivial solution. 
Actually it has been shown that this phenomenon occurs in QED4, if the 
gauge coupling is allowed to be strong. This suggests the existence of the criti- 
cal point /3g such that the nontrivial solution exists only for (3e := 47r^/e^ < (3^ 
(supercritical or strong coupling region). 

The simplest way for obtaining the value of the critical coupling and the 
explicit form of the nontrivial solution is to look for the bifurcation solution of 
the nonlinear integral equation |5^, |56| . The bifurcation solution from the trivial 
one, if any, is obtained from the following |^ linearized equation of Eq. (^3| 



b{x) 



167r2 



d 

dy 



dB{y) 



yB{y) 
y + B\y) 



Kix,y)b{y), 



(9. 



B=0 



namely, 



b{x) 



1 

IGtt^ Jo 



A2 



dyK{x,y)h{y). 



(9.9) 



As the linearized equation does not determine the scale of the solution, we impose 



the following condition to normalize the solution: 6(p^ 



m 



m. The nontrivial 



solution corresponds to the bifurcation solution of the original nonlinear integral 
equation. The critical coupling constant /?g is obtained as the bifurcation point 
of the solution, at which the number of the solution changes. 

Another way of linearizing Eq. (|9.5|) is to study the following equation [H : 



B{x) = mo + 



167r2 



dy 



yBjy) 
y + B^{0) 



K{x,y), 



(9.10) 



where the normalization of the solution is given by m = B{0). 

Both method give the same result at least on the critical point, if they have 
the nontrivial solution: (3^ = 3 for quenched QED4. 



-'^-'^Here we have assumed that the kernel K{x,y) does not include the function B{x). This is 
satisfied at least in the quenched ladder approximation. 
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9.1 strong gauge-coupling region: <C 1 

First, we consider the strong gauge-coupling region, i.e., /? <^ 1. In this region. 



the nonlocal gauge ( |9.2|) can be expanded as follows: 



E 



\ (m + 2)(m + 3) / 



2M2 10 
In accord with this expansion, we have 



(9.ii: 



1 



(9.12) 



Therefore the integrand of the kernel of the integral equation is given by 



_6_ ^ m + 2 



M2 



4 

M2 



6 

~ 8M2 ^ 5 V M2 / 



+ (9.13) 



Making use of the integration formulae 



(»-!)!! 



(n = 0) 



(n+2)!!^ (n = 2,4,6,...) 
(n = 1,3,5,...) 

we can perform the angular integration in the kernel: 



(9.14) 



K{x,y) = AGT 
Then the SD equation reads 



1 - ^/3Gt(x + y) + \(3^Gl{x^ + y^ + ?>xy) 



B{x) =m, + —j^ dy 



yB{y) 
y + B\y) 



K{x,y), 



+ 0{p^). (9.15) 



(9.16) 



where we have defined 



K{x,y) ■.= {AGT)-'K{x,y). 



(9.17) 



Since the kernel K{x, y) is symmetric in x, y, the existence of the nontrivial 
solution is guaranteed by a mathematical theorem for the integral equation with 
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a symmetric kernel ||5^, as long as the cutoff A is finite. However, the existence 
of the solution in the continuum limit A ^ oo is a nontrivial problem which we 
address in this paper. 

In the limit f3 ^ 0, the kernel K{x,y) is a constant: K{x,y) 1. Hence the 
SD equation reduces to the well-known gap equation whose solution is given by 
a constant (i.e., momentum-independent) mass -B(p^) = m: 



m 



mo + 



dy- 



ym 



An^ Jo y + rri 



(9.18) 



In the chiral limit rriQ - 
scaling law: 

1 

9 



0, this gap equation is easily solved to give the following 



47r^ 
GtK^ 



m2 
A2 



In 1 



A2 



This shows that the critical coupling gc is given by 



9c 



(9.19) 



(9.20) 



and that the critical exponent for the dynamical fermion mass has the mean- 
field value, 1/2, apart from the logarithmic correction. To see more details on 
the critical behavior, it is convenient to left mo behind: 



mo , _i , m m^ / A^ 



(9.21) 



This is the equation of state for the Thirring model. The equation of state is able 
to deduce the critical exponents, u, P, 7, a, 6 which take their mean-field values. 
Taking account of the terms up to 0{(3) in the kernel, we obtain 



B(x) 



mo + 



G, 



47r2 



Gt9 



47r2 



9 

-(3Gtx / dy 
8 Jo 



y + B\y) 
yB{y) 



1 



9 



l3GTy 



Hence we have 



B'ix) 



Gt 9 



47r2 



9 /-A' 
^^Gt / dy 
8 Jo 



y + B\y) 
yB{y) 



y + B\y) 



The RHS is an x-independent constant, —ai. This implies that B'{x) = 
as long as B{x) > 0. As B"{x) = 0, the solution is written in the form: 



-aip + m 



m 



ai 2 

—P 

m 



(9.22) 

(9.23) 
-ai < 

(9.24) 
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The solution is determined from the initial conditions: -B(O) = m, i.e., 



m = niQ + 



47r2 Jo y + B'^{y) 



1 - ^PGtV 



and -B'(O) = ai, i.e., 



Gt9 



ai 



47r2 



-(3Gt / dy 



(9.25) 



(9.26) 



Therefore, the coefficients ai and m in Eq. (|9.24| ) are determined by solving 
Eq. (p:^ ) and Eq. (^). 

The critical line in the phase diagram (/?, Gt) is obtained from the bifurcation 
solution h{p^) = —aip^ + m. Hence Eq. (|9.25|) and Eq. (|9.26| ) are respectively 
replaced with 



and 



m = mo + 



ai 



Gt 
4^ 



dyb{y) (l - l(3GTy) 



Gt 9 



47r^ 



9 

/ dyb{y). 
8 JO 



(9.27) 



(9.28) 



When mo = 0, we can obtain only the ratio of two coefficients, ai/m from the 
bifurcation equation, Eq. ( |9.27| ) and Eq. ( |9.28| ). From the Eq. ( |9.28| ), the ratio is 
obtained as 



m 1 + ^/5e^2 ' 

while, in the chiral limit, Eq. ( |9.27| ) leads to 



(9.29) 



m 



1 - jPe9 ' 



where we have introduced the dimensionless coupling constants: 



47r2 GtA^ 



o2 ' 



9 



;,From Eq. 



, we find that the mass function behaves as 

A2 



m 



1 - ^^f^eg'^ 



(9.30) 



(9.31) 



(9.32) 



Note that the coefficients ai and m in the bifurcation solution are different from those of 
the solution of the original nonlinear SD equation, although we use the same notation. 
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Equating two ratios, Eq. ( |9.27| ) and Eq. ( |9.28| ), we can obtain the equation of the 
critical hne. Note that the bifurcation solution gives the exact critical line. 

However this method does not lead to the correct scaling law in the neighbor- 
hood of the critical point of the Thirring model, i.e., in the limit /? — > 0. This 
is similar to the situation encountered in the gauged NJL model: the bifurcation 
solution is insufficient to derive the correct scaling law of the NJL model in the 
e — > limit. Therefore, we adopt the Miransky's prescription to linearize the 
non-linear SD equation to study the scaling law in the strong gauge- coupling 
region /5 ^ 1: 



B(x) = mo + 



G 



T 



47]-2 







y + i?2(0)- 



(9.33) 



By putting x = 0, this SD equation leads to 



mo 



m 



A 



A A3 



yB{y) 
y + 711? 



KiO,y). 



(9.34) 



For the kernel (|9.15| ), it is not difficult to see that the solution is given in the 
form: 



m 



l + 5:(-l)"a, 



2\n 



n=l 



m 



(9.35) 



Once we know all the ratios an/m{n = 1, 2, ...) from the boundary condition, we 
can write the equation of state. 
Making use of the replacement, 



y 



1 



m 



y + m? y + m? 

we can write down the equation of state: 

,3 



mo m m^ 
A A ^ ' A3 



m 
A' 



P. 9 



(9.36) 



(9.37) 



where 



r(/?e,^) ■.= g-^-F{(3,,g):=g-^- T 

JO 



dyB{y) 
A^ m 



K{0,y), 



and 



' dy B{y) 
y + m? m 



K{0,y). 



(9.38) 



(9.39) 
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The general form of A is obtained from 



A 



A2 



dy 



i + E(-ir 



a„A 



2n 



n=l 



y + m 

oo 

i+j2i-iyupGTyy 



m 



y_ 

A2 



1=1 



(9.40) 



It is not difficult to show that A (j^;(3,g^ diverges at most logarithmically for 
large A: 



rrid 



A(—-p,g]r^\n^, 



A 



A2 



m2 



(9.4i: 



since anA?"'/m is finite even in the limit A — > oo. On the other hand, F{(3e,g) is 
finite even in the limit A — > oo. 

In the chiral limit mo 0, the equation of state may have a nontrivial 
solution, that is to say, the fermion mass is dynamically generated when 
T{Pe,g) < (as long as A > 0). Therefore the condition T{Pe,g) = gives 
the critical line in the two-dimensional phase diagram {i3e,g)- Moreover, by the 
standard argument, this equation of state is able to show that the scaling 
law of the gauged Thirring model is given by the mean-field type, i.e., almost all 
the critical exponents defined in the neighborhood of the critical point take their 
mean-field values, apart from the logarithmic correction coming from A. 

If we use the ratio obtained up to 0{Pe) in Eq. ( |9.28| ), we get 



r{Pe,g) = 9 



9 



Hence the critical line is obtained as (see Fig. 6) 
g'^=g-\(ie):=\--Je^O{i?:) = \ 



8 f3^. 



where we have used the critical coupling of pure QED, 

^2 



I.e. 



Air 



IT 

3' 



(9.42) 



(9.43) 



(9.44) 



For example, the dynamical fermion mass obeys the mean-field type scaling: 

m, ^ A[g:\(3,) - g-']'/' (9.45) 

at fixed /3e -C 1- These results should be compared with the result of the effective 
potential in the previous section. The result obtained from the SD equation are 
qualitatively in good agreement with the previous result. 
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It is clear from the above treatment how to include more higher-orders in (3e- 
As all the necessary materials are also presented above, we do not repeat the 
same type of calculation any more. It is not difficult to show that the inclusion 
of higher orders does not change the qualitative feature of the phase structure. 
To make this issue more clear, the detailed numerical studies will be given in a 



subsequent paper ||4^ . 



9.2 strong four-fermion coupling: M <C 1 or Gt ^ 1 

Next, we pay attention to the region of strong four-fermion coupling, Gt ^ 1. 
We introduce a new variable: 



to rewrite the function Dt as 
Then the nonlocal gauge reads 



1 



(9.46) 



(9.47) 



-2 



2^3 



-2 + 



(C-/i 
(C - ^^'f 



J u 



C'-4/i2C + 3/i^ + 2/ln4 



(9.48) 



Note that /x^/C < 1, when k'^ > 0. Expanding the above nonlocal gauge into the 
power-series in /x^/C, we get 



rii^) = 1 - 3^ + 3 f -3 + 2 In 4 
This leads to the following expansion: 



+ 



2\3- 



(9.49) 



+ + 3 3 - 2 In ^ — — + O 



?>D^?{e) + 3M2[d5?)(P)]2 + O [M^[DP{e)f] . 



(9.50) 



Taking into account only the first term in the RHS of this equation, we are lead 
to the SD equation: 



S(p2) = rno + 3e^J 



(27r)4g2 + 52(^2) (p_^)2 + ^2 



(9.51) 
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This is the same SD equation as that treated in the study of the zero-charge 
model . 

First of all, we consider the limit /i ^ with being kept finite and non- 



zero, which corresponds to the g oo limit. In this limit the SD equation ( |9.51 



reduces to the SD equation of the quenched planar QED4 in the Landau gauge: 

B{q^) 1 



Bip") = mo + Se^y 



(27r)4 g2 + 52(^2) (p_^)2- 

The kernel ( ^I6D of the SD equation Eq. (|9.5| ) is given by 

9{x-y) 6{y~x) 



(9.52) 



K{x,y) := — [ (I'd sin^ 

71 Jo [p — qY 



3e^ 



X 



y 



(9.53) 



Note that the nonlocal gauge reduces to the Landau gauge in the limit of quenched 
planar QED4 [Q. The solution of the quenched planar QED4 in the Landau 
gauge is well-known In the chiral limit mo = 0, this equation has a trivial 
solution -B(p^) = for arbitrary coupling. Even in the chiral limit mg = 0, 
however, the SD equation has a nontrivial solution, if the gauge coupling is allowed 
to be strong, i.e., in the region /?e := 47r^/e^ < /3g (supercritical or strong coupling 
region) where the chiral symmetry is spontaneously broken and the fermion mass 
is dynamically generated. 

The simplest way to obtain the value of the critical coupling and the ex- 
plicit form of the nontrivial solution is to solve the bifurcation equation of 
Eq. dH): 



Hp' 



3e' 



(27r)4 g2 ^p_qy 



(9.54) 



The nontrivial solution corresponds to the bifurcation solution of the original 
nonlinear integral equation. The critical coupling constant j3l is obtained as the 
bifurcation point of the solution, /3g = 3. In the neighborhood of the critical point 
I3e = /3g, the nontrivial bifurcation solution in the chiral limit mo = is given by 



m 



m 



1/2 



sm 



1 , f 



+ 5. 



/ sin 5 



pi 



(9.55) 



where 



1, <^p = arctan p. 



(9.56) 



Here the solution is normalized so that h{p^ 



m 



m. 
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This solution leads to the following scaling law of the essential singularity type 
(Miransky scaling) for the dynamical fermion mass m^: 



A 



exp 



rni — 5p 
P 



n 



1,2, 



(9.57) 



The solution with = 1 is the ground state and the solutions with n > 2 
correspond to the excited state in the sense that the n = 1 solution gives the 
lowest (CJT) effective potential than other solutions with n > 2. 

Next, we consider the case of non-zero /i^, but not so large: the ratio 



r := 



A2 



GtA2 



1 



(9.5^ 



is kept finite (and non-zero) in the limit A^ — oo. Therefore this situation exactly 
coincides with that of the zero-charge model treated in [B3. The quenched planar 
QED4 is obtained as a limit — ^ 0. The results are summarized as follows. 
In this case, the kernel (|9.6|) of the SD equation Eq. (|9.5D is given by 



k{x,y) 



77 JO 



3e^ 



6e2 



X + y + + 



1 

1 + - 



X 



P 



y + p 
2 



2^2 



Axy 



\J {x + y + fi'^y — Axy 



(9.59) 



In ref. the approximate analytic solution of the SD equation Eq. ( |9.51 

has been obtained by replacing the kernel with the simplified one: 



K{x,y) = 3e' 



0{x-y) ^ 0{y-x) 



X + fi^ 



y + p^ 



(9.60) 



In the neighborhood of the critical point, the (bifurcation) solution in the chiral 
limit mo = is given by 



m 



»2 + ^2 



sm 



1 £ + p^^. 

2 + fj,'' 



I sin 5^ 



(9.61) 



The equation of state is obtained as follows. 

,2\ 1/2 



mo 



m- 



1 + p Im^^yi.^ 



2p VA^ + yU^ 



sm 



1 1 A!±if!^A' 



(9.62) 
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In the chiral limit mo = 0, we obtain the following scaling law for the dynamical 
fermion mass irid'- 



A 



^1 + r) exp 



P 



- r. 



n 



1,2, 



(9.63) 



It is obvious that, in the limit /x — 0, the solution Eq. ( p.61|) and the scaling 
law Eq. (|9.63|) reduce to those of quenched planar QED4, i.e., Eq. (|9.55|) and 
Eq. ( |9.57|) respectively. 

The critical coupling (or (51) is obtained as the value of parameters at which 
nid/ Is. vanishes. As a result, the critical coupling (or (31) is written as a function 
of r: = e^(r) (or (31 = /3g(r)). For a given r, the critical coupling /3g(r) for 
r < 1 is given by 



:i+p' 



using a solution p of the equation: 



2p 



p2 - 1 



tan 



|ln(l + r-^) 



(9.64) 



(9.65) 



For a given r, if p is a solution of Eq. ( |9.65| ), then so is for —p. Note that p^ > 
for a real value of p. Therefore the real solution of Eq. ( p.65| ) leads to the critical 
point /3g(r) such that /?g(r)//3g(0) < 1. It is shown that the critical value P^i''^) is 
monotonically decreasing in r. 

Note that the critical point should locate on the straight line = r(3e and 
this line coincides with the pure QED axis in the limit r — »• 0, and with the pure 
Thirring- model axis in the limit r 00. Therefore a critical point is obtained 
as an intersection point of two lines, = r(3e and (3e = Pli^) for a given r (see 
Table 1). In Fig. 6, the critical line is drawn by plotting the intersection point 
for various r. Incidentally, the numerical solution of the non-linear SD equation 
Eq. (U) shows [g that 



Mi) 

/5e^(0) 



0.5 



0.17. 



(9.66) 



Moreover, it has been shown that the critical scaling law obeys the mean- 
field type, namely, the critical exponent takes the mean-field value except for the 
possible logarithmic correction: 



md ~ A(e2 - el{T)r\ \og{e^ - e^(r))r', v=- 



(9.67) 
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This imphes that 



ch 



(9.68) 



since the chiral order parameter is related with the mass as follows within the 
approximation Eq. ( |9.6(J| ): 



(9.69) 



The critical exponent for the logarithmic correction in strongly coupled gauge 
theory has not been known up to now to the author's knowledge, although it is 



tried in ref. 56 



It is rather difficult to study the effect of higher orders in analytically 
beyond the lowest order given in this paper. The numerical investigations will be 
given in a subsequent paper 



Table 1 



r 


P / 












1 





0.01 


0.785 


0.619 


0.018 


0.05 


1.020 


0.491 


0.074 


0.10 


1.176 


0.420 


0.126 


0.20 


1.391 


0.341 


0.204 


0.30 


1.557 


0.292 


0.263 


0.40 


1.699 


0.258 


0.309 


0.50 


1.825 


0.231 


0.346 


0.60 


1.941 


0.210 


0.378 


0.70 


2.048 


0.193 


0.404 


0.80 


2.149 


0.178 


0.427 


0.90 


2.244 


0.166 


0.447 


1.00 


2.335 


0.155 


0.465 
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10 Equi-correlation-length line 

10.1 strong gauge-coupling region: /3e <C 1 

Fist we discuss the case of /5e <C 1. /,From the equation of state, Eq. ( |9.37| ), we 
observe that in the chiral hmit mo = 0, ^ satisfies the following equation: 



In other words, the correlation length ^ defined by 

A 



rrid 



10.1) 



(10.2) 



is given as a function of two coupling constants, j3 and g. Therefore the equation 

g-^ - F(/?, g) = -r'A (r^ P, g) , (10.3) 

gives the line of equi-correlation-length C, in the bare parameter space (/5, g). The 
infinite correlation-length limit ^ ^ oo corresponds to the critical line. 



TiP,g) ■.= g-'-FiP,g) = 0, 
since {"M (^"^ P, g) 0. Up to 0(/?), it is obtained by Eq. (|913|) . 



;io.4) 



10.2 strong four-fermion coupling region 

When 1/g <^ 1 and Pe is in the neighborhood of Pe, the equi-^/ := A/m^ line is 



obtained from Eq. (|9.63|) as 



-1 



exp 



2{mT~Sp) 



1 — exp 



2{mT—Sp) 



(10.5) 



This is drawn in Fig. 7 for various values of In the limit ^ — oo, this line gets 
equal to the critical line: 



g ^ = gc\Pe) := Pelexp 



2{nn — 6p) 
P 



1 



(10.6) 



Therefore, in the neighborhood of Pe = P^, the critical line behaves as 

2{n7r-6py 



g;\Pe)=Peexp 



P 



;io.7) 
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All the derivatives of the critical line at (3e = (31 are zero, since the critical point 
(3e = is the essential singularity. 

Here it should be noted that the gauged Thirring model has the dynamically 
generated fermion mass even if 1/(7 < 0. This is interpreted as follows. Al- 
though the original four-fermion coupling is negative for 1/g < 0, namely, the 
four-fermion interaction is repulsive, a pair of fermion and antifermion can be 
condensed to make the bound state due to strong gauge coupling /3e < 1. It is 
known that such a phenomenon occurs also in the gauged NJL model: the dynam- 
ical fermion mass generation is possible even for negative four-fermion coupling, 
if the gauge coupling is sufficiently strong [[15], |T^. 

Therefore this result suggests that the dynamical fermion mass generation is 
possible in the negative four-fermion coupling region 

Gt < G'^iP) < 0, (10.8) 
as well as the positive coupling region: 

Gt > G^tW) > 0. (10.9) 
It is worth recalling the fact ^ that the 1+1 dimensional Thirring model 



indeed exhibits the dynamical mass generation when Gt > and Gt < — vr, i.e., 
G^(0) = and G^(0) = —n. In the gauged Thirring model, the critical values 
G^, G'j. depend on the gauge coupling /3. This issue will deserve further studies. 



10.3 the whole critical line 

We have shown the existence of the critical lines which extend into the interior 
of the phase diagram {P,Gt) from two critical points: (0, Gy) and (/3c, 0). In 
the neighborhood of the respective critical point, it was shown that the scaling 
behavior is consistent with the mean-field theory, except for possible logarithmic 
corrections. We expect that the two critical lines meet together and constitutes 
one critical line as a whole. Then the critical line connects two critical points of 
pure Thirring model and pure QED4. 

In order to show this, we need to solve the SD equation Eq. ( |9.5| ) numerically 
in the nonlocal gauge ( p.2| ) without any further approximation. For this, it is 
necessary to perform double integration over the angular variable and the radial 
variable p^. Such kind of numerical calculations were done in strong coupling 
QED4 beyond the quenched approximation where the angular integration is a 
nontrivial problem, see |^T], The result of numerical calculation will be given 
in a subsequent paper together with other things. 
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11 RG flow and the continuum limit 



11.1 renormalized coupling constant 

By virtue of the reformulation of the Thirring model as a gauge theory, the original 
four-fermion contact interaction has been replaced with the gauge interaction 
mediated by exchanging the gauge boson, see Fig. 8. Therefore, in analogy with 
QED, we can define the renormalized coupling constant between the fermion 
and the gauge boson in the chiral symmetry breaking phase as follows. 



1 1 A;^ 



^ (11.11 



mi 



where the factor 1/47?^ has been introduced to simplify the final expression. Note 
that this definition of the renormalized coupling constant does not depend on the 
longitudinal part of the gauge boson propagator and hence is gauge- independent. 
Q The above equation is rewritten as 

1^ = ?^ ni2) 

11.1.1 strong gauge-coupling limit /3e — 

In the limit /?e 0, Eq. ( |TLID and Eq. (|r|) reduce to 

o^R = 9'§. and '^ = {ang-'f/\ (11.3) 



Substituting this relation into the equation of state Eq. (p.21| ) for the Thirring 



model, we obtain the equation of equi-a^j line (in the presence of the bare fermion 
mass) : 

^ = g{aRg-yl' [g-' - 1 + ' In (l + ang-')] . (11.4) 
In the limit niQ = 0, the an is implicitly given as a function of g and vice versa: 

g-'-l + ang-' In (l + ang-') = 0. (11.5) 

In the chiral symmetry breaking phase, the fermion becomes massive, i.e. the fermion has 
the pole mass irid at = m^. Therefore, we can define the renormalized coupling constant at 
another point: fc^ = 2m^ for fc^ := {p — q)^ = + — 2pqcos-d. However, this changes of the 
definition does not essentially influence the result, see the next subsection. 

We identify the dynamical fermion mass rn^ as a pole mass which should be gauge- 
independent, although we actually use nid obtained from the approximate SD equation in 
the nonlocal gauge. 
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If a value of an is specified, the corresponding value of g and hence md/ is 
uniquely determined according to the equation of state which is derived as a 
consequence of SD equation. In Table 2, the values of g and m^^/A are enumerated 
for various a^. In the phase diagram, an equi-a/j line on which has a specified 
value of an begins at the point on the line /3e = (pure Thirring- model axis), and 
smoothly meets with the piece of the corresponding line in the region f3e ~ 131, 
which is obtained in the next subsection. 



Table 2 





9 


9-' 







1 


1 





0.1 


1.0094 


0.99064 


0.31475 


0.2 


1.0353 


0.96588 


0.43952 


0.3 


1.0739 


0.93118 


0.52854 


0.4 


1.1219 


0.89129 


0.59709 


0.6 


1.2372 


0.80826 


0.69639 


1.0 


1.5085 


0.66289 


0.81418 


2.0 


2.2655 


0.44140 


0.93957 


3.0 


3.0532 


0.32753 


0.99125 



11.1.2 strong four-fermion coupling near j3l 



When (3e is near j3l and 1/g ^ 1, the dynamical fermion mass obeys the 
following scaling law Eq. (|9.63| ): 



m 



^ = (/3e^?)-M(/3e^7 + l)exp 



2{n-n — 5p) 
P 



;ii.6) 



Therefore, substituting Eq. (|11.6| ) into Eq. ( |11.2|) or Eq. (|11.1|) , we get the equa- 
tion for the equi-tt/j line in the region 1/g ^ 1: 



/3e 



OtR^e 



exp 



2{n'K — 6p 
P 



;ii.7) 



In Fig. 9, this line is drawn for various values of a^. The graph shows that 
the line passes through the negative 1/g region for relatively large value of 
{aR > 0.3). Especially, in the limit an 0, the line of equi-au Eq. (|11.7|) reduces 
to Eq. (|10.6|) . As Eq. ( |10.6|) is nothing but the critical line, the renormalized 
coupling constant vanishes on the critical line m^/A = 0. Hence this result seems 
to show that the continuum limit of the cutoff gauged Thirring model is trivial. 
However there is an exceptional point, i.e., {{3^,1/ g) = {P^,0). All the other 
renormalization group flows with non-zero converge to this point. This shows 
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that we can adjust the bare parameters so that the cutoff gauged Thirring model 
moves along the RG ffow with 7^ in the continuum limit and the resulting 
continuum theory obtained in such a way has the nontrivial interaction. 



11.2 constant mass ratio 

In order to get further insight into the renormalized properties of the gauged 
Thirring model, we study the line of equi-mass-ratio between the fermion mass 
and the gauge-boson mass. As the mass of the gauge boson is given by 

m\ = c^Gt^ (11.8) 
within our approximation, the mass ratio reads 

rm:=^ = ^p^ = f3e9^. 11.9 

As the physical mass is gauge-independent, this ratio should be also gauge in- 
dependent. In our treatment, this is approximately satisfied, since we use the 
fermion mass obtained from the approximate solution of the SD equation. 

If we use Eq. (|11.6|) for the dynamical fermion mass, we get the line of equi-r^ 
(mass ratio), see Fig. 10: 

N -1 



g ^ = I3e{{l +rrn) exp 



2(n7r — 6p) 



P 



1 . (11.10) 



The line ( [L1.10| ) with = is equal to the critical line. We will be able to obtain 
the continuum gauged Thirring theory with massive fermion and massive gauge 
boson by approaching the point (/?e, ^/ 9) = {Pe^ 0) along this line with 7^ 0. 
Comparing Eq. (|11.7|) with Eq. ( p.l.l0| ), we can see that the equi-r^ line 



coincides with the eqai-an line, if they are reparameterized according to the 
following relationship: 

1 - PeOiR 

In the limit of approaching the critical point, the above relationship reduces to 
r^n = , ^^^"f , or an = {f3t)-'^^. (11.12) 

Here note that is monotonically increasing in a^. Therefore, for the mass 
ratio to be positive, < < +00, in the continuum limit, the should lie in 
the region: 

0<aR<{P^J-' = ^. (11.13) 
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Therefore, the renormahzed couphng constant an should have an upper bound, 
for the continuum theory to be physically meaningful. ^ In this range of an, the 
equi-ttR lines locate in the positive region, 1/g > 0. 



11.3 scaling of the mass function 

Note that the bifurcation solution ( p.61| ) is rewritten as 



1/2 



where we have used 



sm 



_|_ rp 



1 + r~ 



5n 



I sin 5 



pi 



:ii.i4) 



:ii.i5) 



This shows that the fermion mass function has the universal form on the equi- 
Tra line and hence on the equi-cti? line. Therefore, if we plot the fermion mass 
function h{p'^)/md as a function of p/md-, the curves agree for various different 
values of the cutoff (or rrid/A) on the equi-ttR line. Hence the scaling holds in the 
neighborhood of the critical point (/3e, ^/g) = (Pei^)- 



11 A beyond the quenched approximation 

In this paper we have adopted the quenched approximation to obtain the an- 
alytical solution of the SD equation for the fermion mass function. And, we 
have inevitably substituted the result under our approximations into the m^/A 
of Eq. ( |11.1| ) and ( |11.9| ) to obtain the RG flows. However, the following result 
is unchanged even if we use the exact solution of the SD equation to study the 
RG flows: it is necessary that the four-fermion coupling constant approaches 
infinity g{A) ±oo for obtaining the nontrivial continuum theory in the contin- 
uum limit A — > cx). In fact, if g remains finite in the continuum limit A ^ oo, 
both and go to zero in this limit. Therefore, if we want to obtain the 
nonzero an and r^, g must go to infinity while (3e must converge some finite 
value. Therefore, although the explicit RG flows may change depending on the 
approximation, the above scenario for obtaining the nontrivial continuum limit 
should hold irrespective of the approximation. 
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If we adopt an alternative definition of the renormalized coupling constant, the relation 



(11.11) changes as follows. 

1 



or OR = {(3^y 



However, this does not change any conclusion, e.g. Eq. (11.13) 
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12 Anomalous dimensions 



12.1 strong gauge-coupling region 

First, we obtain the wavefunction renormalization constant defined by 

dniQ 



dm 



12.11 



L 



From the equation of state, Eq. ( |9.37|) , we observe that in the chiral hmit mo = 



r(/3,^?) ■.= g-'-F{P,g) 



A2 V A 

^From the equation of state Eq. ( |9.37| ), we have 



P. 9 



'12.21 



dm. 







m 



dm 



A2 



m 



gT{(3, g) + 3(7— A [-■(3,g\+ g—A [-](3,g 



A 



m 



A3 



m 



A 



;i2.3) 



where the prime in the third term of RHS imphes the differentiation with respect 
to m/A. The third term is calculated as follows. 



A' 



m 
A' 



P,9 



d 



-A 



dm/ A 
-2 



m 



P,9 



m dK^ 
A3 1 



A 



m A^ + m^ m 



K{0,A' 



;i2.4) 



where we have used Eq. 



. Therefore we obtain 



-2gTiP,g)-2g-f 1 + -^ 



m^V' BiA 



'29 



A2 



A2 



'A2 



F{f3,g)-^A 



A2 

rrid 
A 



1 + 



]P,9 



K{0,A'^ 



A2 



K{0,A') 



X 



A{^;P,9 



A2 



-K(0,A^ 



;i2.5) 



In the limit (3 —>■ 0, especially, we obtain for large A/m 



'A2 



1 - 



mi 



A2 



- -1 


, A2 




In^ - 1 







;i2.6) 
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Now we calculate the anomalous dimension of the composite operator ipip 
according to the formula p3l- 



d\nZ„ 



dlnA 

/^From Eq. (|12.5| ), we obtain for large A/m 



A/m^oo 



dlnA 



2 + 



A 



\ A J nid 



n -1 



d 



d\nA 



-ir(0,A2 



'1 



d 



(91nA 



A 



A 



12.7) 



B{A^ 
rrid 



-K{0,A'()12,.8) 



where we have omitted to write the dependence of the coupling constants. Note 
that F{P,g) and ^j^K(0, A'^) are A-independent constants, as discussed in sec- 
tion 9. Hence, the derivative 



, A^ B{A^ 

A^ _|_ ^2 



-K{0,A' 



;i2.9) 



is finite in the limit A ^ oo. This is equivalent to say that A diverges at most 
logarithmically as the cutoff A increases. Therefore the anomalous dimension 7^ 
defined by Eq. ( |12.7| ) approaches the large value in the continuum limit A ^ 00: 



irrespective of the values of the coupling constants, /3,g. 



(12.10) 



12.2 strong four-fermion coupling region 



/^From the equation of state Eq. ( |9.62| ), we obtain the wavefunction renormaliza- 
tion: 



l+p2 



2 + /i^ 



A2 + /i2 ) 



cos 



-pln^- — - + dp 

2 mfi + fx^ 



(12.11: 



Then we have 



' dlnA 



P^~, — 2 



tan 



1 1 A!±i^^A' 



(12.12) 
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If we substitute the scaling law Eq. ( p.63| ) for the dynamical fermion mass 
into the RHS of this equation, the third term vanishes, since the argument of tan 
is n7r(?7, = 1, 2, ...). Hence the anomalous dimension reads 

By putting /i = 0, Eq. ( |12.13| ) correctly reproduces the large anomalous dimension 
of quenched QED4 |F 



7^^ = 1- (12.14) 

If we move the cutoff gauged Thirring model to the critical point in the arbitrary 
way, i.e. irrespective of the RG flows, we get the largest anomalous dimension, 
7m = 2, as in the case of the previous section. However this is physically mean- 
ingless, since this leads to the trivial continuum limit, i.e. ur = 0. 

If the continuum limit A ^ 00 of the gauged Thirring model is taken along 
the RG trajectories obtained in the previous section, the anomalous dimension is 
related to both the renormalized coupling constant and the mass ratio as follows. 

2 + , , 

7v;v = TT = 2-3aR. (12.15) 

This relation implies that the gauged Thirring model has the anomalous dimen- 
sion definitely between 1 and 2 in the continuum limit: 

K < 2- (12.16) 

In other words, the anomalous dimension of the continuum limit of the gauged 
Thirring model is uniquely determined by specifying the RG trajectory along 
which the continuum limit is taken. The equi-7m line is obtained by reparametriz- 
ing the equi-a/j line or equi-r^ line according to Eq. (|12.15|) . The large anomalous 



dimension 7^^ > 1 implies the physical dimension of the four-fermion operator 
gets less than four in 3-1-1 spacetime dimensions: dim[('?/''?/')^] = 2(3 — 7^^) < 4. 
This implies that the four-fermion interaction becomes relevant, which will lead 
to the (nonperturbative) renormalizability of the four-fermion interaction. 
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13 Conclusion and discussion 



In this paper we have proposed a gauge-invariant generahzation of the Thirring 
model, i.e. four-fermion interaction of the current-current type. We call this 
model the gauged Thirring model. The original Thirring model is recovered as 
the strong gauge-coupling limit of the gauged Thirring model. This is in sharp 
contrast with the gauged NJL model in which the NJL model is obtained as the 
weak gauge-coupling limit. 

In the chiral limit mg = 0, it has been shown that there occurs the spon- 
taneous breaking of the chiral symmetry. We have obtained the critical line 
which separates the chiral-symmetry-breaking phase from the chiral-symmetric 
phase. Moreover, we have shown that the phase transition associated with the 
spontaneous chiral-symmetry breaking is of the second order by examining the 
gauge-invariant effective potential for the chiral order parameter. In order to 
study the phase structure in more detail, we have written down the SD equation 
for the fermion propagator in the nonlocal gauge and obtained the bifurcation 
solution of the SD equation in the two regions which are neighborhoods of the 
critical points of the pure Thirring model and the pure QED4. Both approaches 
have given the same values for critical exponents and the same result on the order 
of the chiral transition. 

It is worth recalling that the introduction of the gauge interaction to the NJL 
model of four-fermion interaction has promoted the trivial and perturbatively 
non-renormalizable NJL model into the nontrivial interacting model in the con- 



tinuum limit 0, |12|, [1^, 3. We have searched for such a possibility also in 



the gauged Thirring model. Actually we have obtained the RG flows (lines of 
constant physics) for the gauged Thirring model. Proposing the gauge-invariant 
criterion of triviality or nontriviality, we have examined the continuum limit of the 
cutoff gauged Thirring model and obtained a signal of the nontrivial continuum 
limit. The nontrivial continuum gauged Thirring theory with large anomalous 
dimension between 1 ant 2 will be obtained by approaching the critical point 

g) = {P'^, ±00) along the RG trajectories of constant renormalized coupling. 
This way of approaching the critical point determines uniquely the ratio of the 
fermion mass to the gauge boson mass in the continuum limit. These results are 
based on the approximate solution of the SD equation for the fermion propaga- 
tor. It will be the next subject to study whether the nontriviality of the gauged 
Thirring model survives after the improvement of the approximation. 

There is a possibility of understanding the nontriviality of the gauged Thirring 
theory from another point of view. This is based on the observation that the 
gauged Thirring model may be obtained as a special limit of the perturbatively- 
renormalizable extended model of the gauged Thirring model, in the similar way 
that the gauged NJL model is obtained as a limit of gauge- Higgs- Yukawa system 
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[0. For example, such an extended model will be given by 

+ ^F^''F,^ + ^\{d, - zA,)<j>\' + ml\<j>\' + X{\<j>\r. (13.1) 

The behavior of this model is well known in the weak coupling region: e ^ ^ 
1) in which the perturbation theory in the coupling constant is applicable well, 
see [^, However, one does not know much of the behavior of this model 
in the strong coupling region: e ^ 1{(3 <^ 1). For the gauged Thirring model, 
therefore, the usual RG function obtained up to some finite order of perturbation 
theory is not necessarily powerful in sharp contrast with the gauged NJL model, 
since we are dealing with the very strong gauge-coupling region. Therefore, we 
need some non-perturbative methods. 

The lattice theory is one of the possibilities. In relation with this, we must 
mention the recent work of Jersak et al. pO| in which a gauge-invariant gener- 
alization of the NJL model was proposed. In our opinion, their model is nothing 
but the lattice version of the gauged Thirring model. They start from the lattice 
formulation with the compact gauge field which will be necessary to extend the 
model into the non-Abelian case. However, their result shows very complicated 
phase structures which make the analysis of the continuum limit quite difficult. 
This may be an artifact coming from the compact formulation of the Abelian 
gauge theory and hence the comparison with the non-compact formulation will 
be needed on the lattice. However, their studies are quite interesting, if their 
model is identified with a simplified version of the non-Abelian gauged Thirring 
model. 

Actually the extension of the gauged Thirring model to the non-Abelian case 
is quite important, since such a model will be a low energy effective theory of 
QCD. Our formulation of the gauged Thirring model is superior to the previous 
formulation ^] in a point that our formulation respects the gauge invariance 
from the beginning as well as all the symmetries of the original QCD. The analysis 
of the non-Abelian case will be given in one of the subsequent works 



Another purpose of this work is to clarify the gauge-invariant mechanism for 
mass generation so that it can be replaced with the usual Higgs mechanism. In 
the very massive case mo ^ 1, we can shed light on this problem through the 
procedure of bosonization. Actually, the 2-1-1 dimensional massive Thirring model 
can be bosonized into the Maxwell-Chern-Simons theory in which the gauge boson 
mass is generated through the Chern-Simons term without breaking the gauge 
invariance [^4|, This scenario can be extended into the higher dimensional 



case In 3-1-1 dimensional case, it is shown |27| that the Thirring model 



can be mapped into the antisymmetric tensor gauge theory. And the mass for 
the tensor gauge field is given in the gauge-independent way. However, we need 



56 



further studies to clarify the non-perturbative renormahzabihty of the tensor 
gauge theory preserving the unitary. This subject will be discussed elsewhere. 

Finally, we mention other possibilities of order parameters. It is very inter- 
esting to study the gauge-invariant (neutral) fermion field F{x) := (f)^{x)ilj{x) 
and F{x) := ■?/^(x)0(a;). Although the neutral fermion gives the same local order 
parameter as the charged fermion ip: 

{F{x)F{x)) = {ij{x)^{x)), (13.2) 

they may have the different two-point functions: 

{F{x)Fiy)) = {i,{x)<P{x)<p\y)ij{y)). (13.3) 

This suggests that the neutral fermion mass mp niay be used as a different non- 
local order parameter. If we regard that mp = implies = 0, we can see that 
7^ implies mp ^ 0. However, = leads to two possibilities: rriF = 0, 
or rriF ^ 0. Hence the chiral-symmetric phase with {ilj{x)ip{x)) = can be dis- 



tinguished by the neutral fermion mass ttif, as suggested by Jersak et al. p9 



Moreover, it is also interesting to make clear the interweaving of the chiral and 



Higgs phase transition There are a lot of questions which are waiting to be 



clarified in the future works 



Acknowledgments 

The author would like to thank Gerrit Schierholz and all members of HLRZ for 
warm hospitality in HLRZ, Forschungzentrum KFA Jiilich as a research fellow 
of the Alexander von Humboldt Foundation. He is also very grateful to Jiri 
Jersak for valuable discussions on the dynamical fermion mass generation by 
strong gauge interaction shielded by a scalar field. This work is supported in part 
by the Alexander von Humboldt Foundation and the Grant-in-Aid for Scientific 
Research from the Ministry of Education, Science and Culture (No. 07640377). 



57 



A Inversion method 
A.l general strategy 

First we add the source term Sj into the original action So := J d^xCo[(j)i{x)] 
where 0i denotes various fields collectively. For example, in order to calculate the 
local order parameter (fij{x) :— {0{x)), we choose the source term: 

Sj := J d^xCj = j d^xJ{x)0{x), (A.l) 

where J is the artificial source. Then we calculate the vacuum action (generating 
functional of the connected Green's functions): 

W[J] = - In Z[J]. (A.2) 

where Z[J] is the partition function in the presence of the source: 

Z[J] := J [dcj>i{x)] exp [-{So[cj>i] + Sj[0])] . (A.3) 

The order parameter 'fj{x) is calculated for the theory with the action: S — 
So + Sj: 

^^^^^ " JjU- ^^-'^ 

It should be noted that the order parameter ipj (as a function of J) is non-zero 
for J 7^ and vanishes in the limit J — > at this stage. 

The effective action T[ip] is defined through the Legendre transformation: 

r[ip] := W[J] - J d^xJ{x)ipj{x). (A.5) 

Thus we have the relation: 

The artificially introduced source J has to vanish in order to recover the original 
theory, i.e. 

0=-|^M. (A.7) 

S(fi{x) ^ ' 

If the operator O has been chosen so as to break the symmetry of JC, the nontrivial 
solution (/? 7^ of this equation gives the spontaneous symmetry breaking solution. 
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For the translation-invariant order parameter, we can define the effective potential 
V{v) by 

a..^l,o^. (A.8) 

Note that the local order parameter is allowed to be a composite operator, e.g. 

0{x) = ^(x)^(x). (A.9) 
This case has been treated in sections 4 and 5. 

A. 2 inversion formula 

Suppose that the order parameter ip is a. function of the source (i.e. probe) J and 
is calculated in a perturbative power series in g: 

oo 

^ = m = j:9''fn[j]. (A.io) 

n=0 

This is called the original series. We know how to calculate each term of the RHS 
of Eq. ( |A.1CI| ), since this is nothing but the perturbation series. 

We need the inverted function of Eq. ( |A.10|) , i.e. J is a function of ip: 

j = r'H (A.ii) 

We frequently encounter the situation in which it impossible to obtain the exact 
form of the inverse function explicitly. However, we can get the above function 
( |A.11|) by a perturbative power series of g: 



n=0 



which we call the inverted series. The important point is that in Eq. ( [A.12[ ) ip 



is regarded as the order of unity. Substituting Eq. (|A.12|) into Eq. ( |A.10| ) and 
expanding the RHS in a power series in g, we obtain 

oo oo 
71=0 m=0 

= fo[ho[^]+ 9hi[(p]+ g'^h2[(p] + ...] 

+gfi[ho[^] + gh[^] + ...]+ //2[/ioM + ghl^] + ...] + ... 
= fo[hoM + g{fl,Mh[^] + MhoM} 

+9'{fo[hoMh2[^] + \UhoMhM 

+f[[h,[^]]h,[^] + f,[h,[^]]} + ..., (A.13) 
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where the prime denotes the differentiation with respect to J . On the RHS of 
Eq. ( |A.13| ), if is regarded as order unity so that ip on the LHS of Eq. (|A.13|) 
should be regarded as order unity. Then we can determine hm[^] by equating 
each power of g on both sides of Eq. (|A.13|) and get the inversion formula: 



hi[ip\ 
h2[p\ 



'm 
- 1 fit 



J=ho[^p\ 



f^[J]hl[v]+f[[J]h,[p]+f,[J] 



(A.14) 



J=ho[ip\ 



Once we know the inverted series Eq. (|A.12|) according to the inversion formula 
Eq. (|A.14 ), the effective action is obtained from Eq. ( |A.6| ) as follows. 



d(f)h[(f)\ 



n=0 



d(j)hn[4>\. 



(A.15) 



A. 3 nonlocal source 



The inversion method can be applied to the nonlocal order parameter as well as 
the local one. For example, in order to calculate the nonlocal order parameter 
ip{x,y) := {0{x,y)), We can add the source action 



Sj = j d^x J d^yj{x,y)0{x,y). 

In this case, we have 

EM := W[J] - fd^'x f d''yJ{x,yMx,y), 



(A.16) 



(A.17) 



and 



¥^{x,y) = 



J{x,y) 



SW[J] 
6J{x,y)' 



6ip{x,y)' 

If we choose the fermion propagator as a nonlocal order parameter: 

ip{x,y) := {il){x)i){y)) := S{x,y), 



(A.18) 



(A.19) 



(A.20) 
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the lowest order inversion in the couphng constant gives the Schwinger-Dyson 
equation in the ladder approximation. Actually, in QED, we get (see [^) in 
Minkowski spacetime: 

J{p) = tS-\p) - tS^\p) - ze' I ^rS{p + k)YD^S}{k), (A.21) 

where D^^}{k) is the bare photon propagator. If we set J = 0, this coincides 

with the ladder SD equation for the fermion propagator. Using the following 
identification: 

m - -|||. (A.22) 



we get the CJT effective action 

T[S] = -^ J -^tr[\n S{p)-S{p)S^' 



(2vr)^ 



This is the case treated in sections 7 and 8. 
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Figure Captions 



Figure 1: The schematic phase diagram of the gauged Thirring model in 3+1 dimen- 
sions at mo = 0. 

Figure 2: The shape of the effective potential when r — 0. (a) Thirring4 model, (b) 
QED4, (c) gauged Thirring4 model. 

Figure 3: The critical line in the phase diagram k) obtained from the gauge- 
invariant effective potential. 

Figure 4: The graphical representation for (a) the CJT effective action, and (b) the 
SD equation obtained from the CJT effective action. 

Figure 5: Plot of the nonlocal gauge-fixing function r]{k'^). 

Figure 6: The critical line in the phase diagram (/?e, obtained from the Schwinger- 
Dyson equation. 

Figure 7: Equi-.^ (correlation length) or equi-md/A (fermion dynamical mass) line in 
the phase diagram (/3e, 1/ g) obtained from the Schwinger- Dyson equation. 
The line is parameterized from above to below as :— m^/A = 0.0, 0.05, 
0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.40. 

Figure 8: Definition of the renormalized coupling constant aR. 

Figure 9: The equi-cti^ (renormalized coupling) line in the phase diagram [(5f.,l/g) 
obtained from the Schwinger- Dyson equation. The line is parameterized 
from above to below as ccr = 0.0, 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 1.0, 1.5, 2.0, 
2.5, 3.0. 

Figure 10: The equi-r^ (mass ratio) line in the phase diagram (/?,., 1/g) obtained from 
the Schwinger-Dyson equation. The line is parameterized from above to 
below as = 0.0, 0.5, 1, 2, 3, 10, 100. 
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